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MOLECULAR VIBRATIONS AND MOLECULAR DIMENSIONS IN 
SOME LARGE MOLECULAR SYSTEMS 

B.L. DEOPURA 
Ph.D. 

Department of Physics 
Indian Institute of Technology Kanpur 
June 1 97 1 

The thesis is concerned with the vibrational studies 
of 1 , 3 j5-Triamino -2 ,4 ,6-Trinitrobenzene (TATNB) and a study 
of low-angle X-ray scattering from polyethylene before and 
after irradiation. In the former a con^lete normal mode 
analysis including evaluation of thermodynamic quantities 
and mean square amplitudes has been carried out. Infrared 
absorption, Raman scattering and inelastic neutron scattering 
are used as tools of molecular spectroscopy. The X-ray 
scattering at low-angles from polyethylene is analysed on the 
basis of a two-phase model and position correlation functions 
are obtained. The peak positions are con^ared with the 
electron microscopic observations. 

The first chapter forms a general introduction to the 
thesis, the importance of vibrational analysis in regard to 
the assignments of the observed frequencies and detailed 
spatial nature of the vibrations of a molecule is discussed. 
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Such studies are very useful in case of macro molecular systems 
•with various conformations. The study of low-angle scattering is 
important for t-wo main reasons. First, it can provide 
information on how the macromoleculax chains are organised 
into the large scale units. Second, in some cases it may be 
possible to obtain information about the internal structixre. 
Further , lo-w - angle X-ray scattering could give information 
regarding the crystallite size, orientation, degree of order 
etc . 

The second chapter of the thesis is concerned with 
the theoretical approach used for the treatment of the normal 
vibrations of a finite molecule^ It is based on Wilson' s GF 
matrix method. The usefulness of group theoretical principles 
in factoring the dynamical matrix and in classification of 
various normal modes is discussed. In addition to infrared 
absorption and Raman scattering, inelastic neutron scattering 
is also an inportant tool of molecular spectroscopy. It is 
because inelastic neutron scattering is not limited by symmetry 
dependent selection rules and as such in principle should 
provide 'complete' information about a vibrating system. The 
theoretical principles involve<l are outlined in this chapter - 

The third chapter deals with the analysis of the TATNB 
molecule based on Wilson's method. Factor group analysis is 
based on point group symmetry. The assignments of the 

observed frequencies are made on the basis of potential energy 
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distribution, symmetry principles, peak intensities, bond 
orders and g?roup frequencies. In neutron spectrum the fact 
that the intensities of various peaks are wei^ted by the 
square of the amplitude of motion of the protons is used. 

NH^ and NO^ group modes, X-sensitive modes and ring modes 
are discussed in details. Urey-Bradley force field is 
obtained. Spatial representation of a few important normal 
modes is presented. The thermodynamic parameters are evaluated. 

The next chapter gives the theoretical details regarding 
X-ray scattering studies at low-angles. A two-phase electron 
density model for semicrystalline polyethylene is presented, 

A layer-like structure is considered. Theoretical correlation 
function from the autoconvolution of the electron density 
v^ariations in the san^le-model is calculated. It can be 
compared with the experimental correlation function obtained 
from fourier transform of the scattering curve. The variables 
of the model (crystalline and amorphous lengths and the widths 
of the distribution functions) are varied to match the 
theoretical curve with the experimental one. 

In the fifth chapter experimental and theoretical studies 
on lew - angle scattering from a polyethylene sample which has 
been annealed just below the melting point and cooled very 
slowly to room ten^Jerature are presented. Luzzati-Baro camera 
is used and for slit hei^t correction, the slit was regarded 
as infinite. The position of the first taa.x±xaxm at 118 a.u. 
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agrees with the electron microscopic observations. In order 
to test the validity of a clear-cut two phase model, studies 
are made on neutron irradiated sample . Since irradiation 
brings about cross-linking, the analysis of the experimental 
curve is complicated because the cross-linked regions could 
have a higher electron density and the distribution of layer 
thicknesses in two phases is broadened. Finally, theoretical 
calculations for a layer -structure and two phase model are mad© 
to fit the experimental correlation curves. Normal distribution 
functions for the average lengths of the crystalline and 
amorphous regions are used. The model parameters are reported. 
Preliminary observations on low-angle X-ray scattering from 
wool are also given. 

The last chapter of the thesis deals with the approxi- 
mations involved in the work and outlines a programme for 
further work. Since TATNB is insoluble in almost all the 
'o'* vents, the spectra are obtained in the solid state. This 
necessitates calculations for a three dimensional system 
taking into account all intermol ecu lar interactions and further 
one has to consider space group symmetry fci' determining optical 
activity and symmetry speciws of different moaes. For low- 
angle scattering the presenr two phase model is certainly not 
the best because in practice a clear-cut demarcation between 
ordered and disordered phase is impossible to make. It is 
expected that a model with more phases would yield better results 
These and several other limitations are discussed in this 
concluding part of the thesis. 
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CHAPTER - I 
INTRODUCTION 

A study of the normal modes of vibrations of a 
molecule is interesting not only in knowing the nature of 
the vibrations taking place but also in understanding the 
thermodynamic properties of the system. The numerical value 
of a vibrational frequency is chiefly a function of the 
magnitude of binding (force constants) among the constituent 
atoms and their masses, whereas the physical nature of a 
vibration depends mainly on the symmetry properties of the 
molecule. If the symmetry and roughly the nature of bonding, 
not the magnitude,, are known in a molecule, the physical nature 
of the normal vibrations can be understood, and thus some of 
the spatial characteristics can be elucidated, like the 
direction in which electric moment, if any, changes during 
a vibration and therefore the selection rules, the degeneracy 
of the normal modes etc. 

Infrared absorption and Raman scattering are mostly 
used as tools of molecular spectroscopy and the assignments 
of the observed peaks without normal mode calculations, is 
done primarily on the basis of group frequencies and thus 
remains somewhat arbitrary. It may not be so bad for small 
molecules. The experimental information itself is generally 
incomplete because several peaks may not be observed either 
because symmetry dependent selection rules do not allow 
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tiiem or else they are too weak to be seen in the spectrum. 

A complete vibrational analysis of a molecule leads to all 
the internal modes and it is possible to obtain details 
about the nature of displacements of various atoms in each 
mode^’^. More important, tlie potential energy distribution 
in various degrees of freedom enables a definite assignment 
of the modes observed. The assignments of the group frequencies 
in model-compounds and identical environments are also very 
helpful in such studies. 

Recently, much work has been done on the vibrational 

3-5 

analysis of long-chain helical macr omolecules , where the 

vibrations of adjacent units are related through a phase- 
factor. Conformational sensitive modes and their dispersion 
are calculated. From the dispersion curves a one dimensional 
phonon distribution function is obtained. 

For 1 , 3 , 5 -triamino -2 ,4,6-trinitroben2;ene (TATNB) , 

molecules such as s-trini trobenzene , aniline, p-nitr oaniline 

could be regarded as model-compounds. They provide useful 

initial information for vibrational analysis of TATNB. There 

are, however, several ambiguities and inconsistencies in the 

assignments of peak positions observed in infrared absorption 

6 —8 

and Raman scattering in these systems . These and trans** 
ferability of force constants are discussed in details in 


Chapter III. 
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1 1 

A study of small-angle X-ray scattering' from 
macromolecular systems is of importance for two main reasons. 
First, it provides information on the organization of macro- 
molecular chains and secondly, in some cases it is possible 
to obtain information about the internal structure of the 
element associated with the small-angle scattering, since 
the transform of the element may determine some features of 
this scattering. For exan^sle , the discrete equatorial 
diffraction maxima can arise from either inter -unit inter- 
ference or intra-unit interference in a unit containing a 
small number of fibrous elements. In the first case the 
distance between the centres of units is assumed large enough 
that the first order diffraction raaxi:rjum occurs at a value of 
(2 Sin &/ A) , which is of the order of l/r where r is the 
distance between the centres of the units, 2© is the scatter- 
ing angle and A is the wave length of the radiation. In 
the second case units are considered to scatter independently 
and the maxima at small angles arise from the lack of extended 
lateral order. 

rfith the common small-angle scattering apparatus, 
the scattered intensity can be measured only on a relative 
scale. A great improvement is achieved when the scattered 
beam is measured on an absolute scale i.e. when the scattered 
intensity is compared to the energy of the primary beam. In 
this case it is possible to determine several additional 
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parameters, such as molecular -weight and internal hydration 

in the case of a solution of globular partxcles * For 

helical macromolecules in solution which behave 3-ike rigid 

rods it is possible to determine linear mass of the rods and 

repeat distance^ Because of the great interest in linear 

polypeptides and nuclic-acid structure, such measurements 

are most fruitful as regard to their conformational analysis^ 

Small-angle X-ray scattering studies cannot provide 

any detailed and exhaustive description of the system and 

one has to contend with average characteristic parameters# 

Specially, in cases where polydisper sity or high anisotropy 

of the particle is present, one can obtain only a statistical 

average for the size and shape* 

In the work outlined in Chapter IV and V, small-angle 

X-ray scattering from a well annealed polyethylene sample 

is presented# Small-angle studies for polyethylene are 

reported by several earlier workers # However , their 

work has been mostly confined to measurement of long periods 

1 9 20 

in fibrous and melt-crystallized samples. Recently, Vonk ^ 
et al# have reported a study of correlation functions based 
on a two -phase parallel layer model with different average 
electron densities# However, their results vary within wide 
limits# The problem has been restudied by carefully taking 
into account the history of the sample# The experimental 
correlation function is obtained by fourier transforming the 



scattering curve. The position correlation functions have 
an important advantage over the direct analysis of the 
scattering cixrve . Interpar tide interferences are auto- 
matically taken care of in the former approach, 
effects are very important in solid state. 


These 
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CHAPTER - II 


THEORY OF MOLECULAR VIBRATIONS 

MOLECULAR VIBRATIONS 

In this section a treatment of the molecular dynamics 
of finite molecules is presented. A molecule with N atoms 
has (3N-6) or (3N-5) vibrational degrees of freedom depending 
on whether it is non-linear or linear. The six non-genuine 
degrees of freedom are associated with free translations and 
rotations of the molecule. For a linear molecule, the rota- 
tional degrees of freedom are only two. Corresponding to 
these (3N-6) or {3N-5) vibrational degrees of freedom, there 
are as many fundamental normal modes of vibration. All other 
modes are linear combinations of them. The Wilson GF matrix 
method is mostly used for calculating the normal modes of 
vibration. In this method ’internal coordinates' are used to 
describe the motion; these are changes in bond lengths, bond 
angles; out-of -plane wags and torsional motions. Physical 
nature of the vibration can be easily visualized in terms 
of internal coordinates. It also enables one to transfer 
force constants from a group in a molecule to a group in 
another in similar environment . The molecular vibrations 
are assumed to be small and harmonic and thus transformations 
between various displacement coordinates are linear. They 
are easily expressed in matrix notation. The transformation 
of cartesian coordinates X to internal coordinates R is 
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expressed as 


R = BX 


( 2 . 1 ) 


B is the matrix of transformation. The kinetic 
energy in terms of internal coordinates is given by 

2T = R' G“^ R (2.2) 

R' is the transpose of R. The G — matrix or inverse 
kinetic energy matrix in terms of B matrix is written as 

®ki ®li / “i = 1,2,...,3N-6) 

i=1 .. (2.3) 

where m^ is the mass of the i atom, n is the number 
of atoms. 

Potential energy is given by the relation 



2V = r' FR 


(2.4) 


where F is the potential energy matrix and is related to 
force constants through Z-matrix by the relation^ 


F 


ij 



Z. f, 
X jk k 


( 2 . 5 ) 


The problem of molecular vibrations then reduces 
to the following secular equation 


J GHF - EA I = O 


( 2 . 6 ) 



•*0 


J \' s are the eigenvalues relate to the vibrational 
frequencies V hy the relation = 4 n ^ c . Normal 

iC ii* 

coordinates Q are related to internal coordinates R by 
the matrix L as 


R = LQ (2.7) 

where the column vectors L are given by the relation 

(GF)L = XL (2.8) 

The solution of the above equation yields the eigen- 
vector L of the dynamical matrix, 

A measure of the nature of the normal vibrations 
(mixing of various modes) is the fractional potential energy 
associated with each internal coordirate. The fractional 
potential energy of the k normal mode associated with 
i^^ internal coordinate is given by^ 

(P.E.)^ = L., F.. /X, (2.9) 

The potential energy distribution is useful in 
making assignments of various frequencies. 

SYMMETRY PRINCIPLES 

The in^ortance of the symmetry of the molecule and 
group theoretical ideas may be mentioned at this point , The 
symmetry and geometry of the molecular -model can be used to 
determine the number of normal modes belonging to different 



symmetry species of the molecular point group, selection rules 
for infrared absorption and Raman scattering, factorization 
of the secular equation and other useful information. 


The number of normal modes in a given symmetry species 

2 


of the point group is given by the relation' 


OO* 

-R -^R 


( 2 . 10 ) 


over 


.y, 



g is the number of symmetry operations and sum is 

all such operations. is the character of the 

K. 

symmetry species, for the operatfon R» Characters 
can be written as 




R 




( 2 . 11 ) 


in which B-re the diagonal coefficient s, in the 

transformation 


t 




( 2 . 12 ) 


which corresponds to symmetry operation R. q. and q. 

are the displacement coordinates before and after the 

symmetry operation. From relation (2.11) and (2,12) it 

follows that can be found by the number of unshifted 

2 

atoms with an appropriate multiplicative factor of the 
symmetry operation. 

Symmetry considerations can also predict the modes 
which will be allowed in infrared absorption and Raman 
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scattering. In the former case, the intensity is governed 
by the changes in dipole moment whereas in the latter case, 
it is dependent on changes in polarizability. In a given 
symmetry point group the three components of the dipole 
moment transforms according to one of the species and the 
normal modes belonging to these species will be infrared active 
with the corresponding polarizations. Similar consideration 
applies for changes in the six polarizability components and 
Raman activity of the normal modes. 

As mentioned earlier, the symmetry of the molecule 
can be exploited in solving the secular equation by trans- 
forming the kinetic and potential energy matrices in the 
space of symmetry coordinates. These are written as linear 
combinations of the internal coordinates with coefficients 
obtained from the Character Table of the molecular point 
group. For a non-degenerate symmetry species the symmetry 
coordinate can be written down as 

s(y) ^ (2.13) 

R 

is same as in equation (2.10) . Notation RS^ stands 
R 

for the coordinate, to which the displacement of the origi- 
nally used coordinate is transformed by the operation 

R. Sum is over all the symmetry operations. N is the 
normalization factor. Each symmetry coordinate belongs to 
one of the symmetry species and in a way represents a normal 
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mode. The dynamical matrix when transformed in the space of 
symmetry coordinates, factors in a block diagonal structure, 
each of these blocks belongs to one of the symmetry species. 

I f 

Let U represent the transformation matrix between internal 
coordinates R and symmetry coordinates S then 

S = UR (2.1^) 


The U-matrix is used po get the block diagonal form 
of G and F matrices by similarity transformation; 


G® = U' G U 
F® = U' F U 


(2.15) 


U is a real matrix, u' is the transpose of U. 
The secular equation (2.6) reduces to the form 
|Q.(X)j'()^)_£7Vj _ 0. These can be solved completely 

independent of each other for individual symmetry species. 


MEAN SQUARE AMPLITUDES 

A mean square amplitude matrix ^ is defined as 


■£ = R R"^ (2.l6) 

R is a set of coordinates for the molecule . R^ 
represents the trans-conjugate of R. In general, a mean 
square an^litude matrix (2.l6) can be constructed using the 
coordinates which need not necessarily be complete or 
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independent . A set of conventional valence coordinates can 
be fruitfully used. The transformation matrix L between 
normal coordinates Q and internal coordinates R can be 
expressed in the form: 


^i " ^ik \ 

k 


(2.17) 


The ^ matrix equation (2.l6) expressed in terms 


of normal coordinates is 



1= 

(2.18) 

or 

1= LbL- 

(2.19) 

where 

fi = > 



using the harmonic oscillator model, we get 

Q > = — Goth (hc<v'/2kT) (2.20) 

4 TT'^c io 

w is the frequency of the oscillator, T is the absolute 
ten^jerature , The elements of ^ matrix for a given normal 
mode k can be written down as 



Y" * 

/ ^ik ^k 

k 


( 2 . 21 ) 


These quantities are independent of the special set 


of internal coordinates which might have been used in the 
normal coordinate analysis for the actual determination of 
L coefficients. At absolute zero temperature 
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Sv(0) 


rr c ui. 


(2.22) 


r (0) 

i j 


L., L. . 
xk i j 


(2.23) 


The corresponding mean-square amplitudes of vibration 
■will be found on the main diagonal of the 5^ matrix. 

ifMm- 

THERMODYNAmC PARAMETERS 

For calculating the thermodynamic parameters the 
rooleeules are assumed to be free. The rigid-rotor and harmonic 
approximations are assumed to be valid for rotations and 
vibrations respectively. The expressions used for calcula- 
tion of thermodynamic quantities^ are as follows; 


= 4r + (g) R ) i i == 

^ ! 1 - exp(- w^ hc/kT)j 


d. w. exp(- w. hc/kT) 

' F ~-H2 (2*24) 


Heat Content 


H°-E° 

o 


4r h- R(^) — 


^i ^i w^^c/kT) 


exp (- w. hc/kT) 


(2.25) 


S° = Entropy (v) + Entropy (Tr) + Entropy (Ro) (2,26) 
Entropy (v) = R ^ log^^1-exp(- w^hc/kT)^ + 


d^ w^ exp (-w^ hc/kT) 
^ 1 - exp(-w. hc/kT) 


(2,27) 


X 
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Entropy (Tr) 


r 

I 

i ■ 


2.2868 I 5 log^Q(T) + 3 log^^CM) j -2.3135 

(2.28) 

Entropy (Ro) = 2.2868 | 3 log^^(T) + 3 ^og^^iA B C) 


- 2 log^ (iT") j - 5*3838 


( 2 . 29 ) 


Free Energy 
F°-E° 


Entropy(Tr) + Entropy (Ro) 


z 


log. 


1-exp(-w^ hc/kT) j - 7*9^52 ( 2 . 30 ) 


Summations are over all the vibrational frequencies, 
is the degeneracy of the vibrational mode. E° is the 
zero point energy. The notations used in above expressions 
are given below; 

T 


w . = 

1 


Temperature in Absolute scale 
till 

i Vibrational frequency 

. th 


d^ = Degeneracy of the i mode 

¥ = Molecular weight 

^ = Symmetry Number ( .j- = 6 for point group) 

A, B, C = Moment of inertia in atomic mass units and 
angstrom units . 

Entropy(v) = Entropy of the molecule due to vibrations 

Entropy (Tr) = Entropy of the molecule due to Translation 

Entropy (Ro) = Entropy of the molecule due to Rotations 

R = Gas constant 

h = Plank's constant 

c = Velocity of light. 



THE USE OF NEUTRON SCATTERING 


The technique of inelastic neutron scattering has 

been applied to study the low lying vibrations in molecular 
7-10 

systems • The usefulness of this technique arises from 
the fact that the energy 'of the thermal neutron is of the 
order of low frequency vibrational and lattice modes and 
their wavelength is comparable to tha^t of the atomic spacing. 
This results in energy and momentum transfer of the thermal 
neutrons to be readily observed. Since neutron scattering 
involves short-range neutron-nucleus interactions in contrast 
to the electromagnetic interactions in infrared absorption 
and Raman scattering, the resulting intensity does not 
depend on the changes in electric dipole moment and polari- 
zability, Because of the symmetry dependent selection rules, 
only a part of the information regarding molecular dynamics 
can be obtained from infrared absorption and Raman scattering, 
whereas neutron scattering being free from such constraints 
provides additional information about molecular systems. 
Experimental Technique 

The neutron scattering experiments have been obtained 

from a time -of -flight spectrometer details of which are 

8 

available elsewhere . A short description of the spectro- 
meter is given here. A qua si -monochromatic beam of neutrons 
with a cut-off in thermal energy region is obtained by 
passing through a polycrystalline Beijrilium sample. Neutrons 
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scattered by the sample S in the direction at right angles 
with respect to incident beam are chopped into brusts by a 
rotating collimator. The chopper thus provides a reference 
time at which the neutrons initiate their flight to the 
detector. Their velocity and hence energy distribution is 
determined by measuring their flight time over a 5 -meter 
path between chopper and detector. Figui’e 1 shows the sche- 
matic drawing of the time-of -flight spec trometer - 
Theory 

The scattering of neutrons can be coherent or inco- 
herent phenomena depending upon the orientation of nuclear 
spin and isotopic effects. The scattering of neutron wave 
from a system of identical nuclei of zero spin would be, 
except for magnetic int orations with atomic electrons, entirely 
coherent and interference among the scattered waves is possible. 
For a system of identical nuclei with nonzero spin, the neutron- 
nuclear interaction depends on the relative orientation of 
neutrons and nuclear spins which gives rise to incoherence. 

In this case, each nucleus scatters independently. For nuclei 
with zero spin, incoherence can persist if the scattering 
length fluctuates from nucleus to nucleus because of isotopic 
effects. Coherent and incoherent scattering together allow 
a complete investigation of the structure and dynamics of the 
scatter er. For example, elastic coherent neutron scattering 
is used for structural investigations and inelastic coherent 
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neutron scattering can be used f'or the deternination of 
dispersion relations for oriented single crystals. On the 
other hand, incoherent (elastic and inelastic) neutron 
scattering can be used for studying the dynamics of the 
atomic system. Several molecular systems contain appreciable 
amount of hydrogen atoms which gives rise to a very large 
incoherent scattering cross section for neutrons. Considerable 
information for molecular vibrations of such systems can be 
derived. Since the cross section for hydrogen is at least 
an order of magnitude greater than that for other elements, 
the scattering from such systems is largely incoherent , We 
can assume that neutron scattering in hydrogenous systems 
is essentially due to hydrogens and incoherent in nature , 

¥e will, therefore, limit our attention to the inelastic 
incoherent scattering of neutrons. 

The fact that the interaction between the neutrons 
and the phonons is of short-range and weak, allows the use 
of the Fermi psoudo -potential. 

U (f , Rjj) = ^ a^ &{f - Rj^) (2.31) 

1 = 1 

where m is the neutron mass, r and R are neutron and 
nuclear position vectors and the scatterer is a system of 
N nuclei each labelled by subscript 'l' . The Fermi 
pseudo-potential contains as a parameter the bound nucleus 
scattering length 'a' . This is a measure of the strength 
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of the interaction. For a scattering: system of identical 
nuclei having zero spin, ’a’ is simply a constant. If the 
scatterer is an isotopic mixture or if the nuclei have non- 
zero spin, the scattering length has to be modified to talce 
into account the spin dependence of the interaction- The 
effect arises because of the interactions in which neutron 
and nuclear spins are parallel, are different from those in 
which the spins are antiparallel. The total interaction in 
terms of the neutron and nuclear spin operators is then 


taken to be 

.. (2.32) 

where b^ is the spin dependent scattering length 1 and 
are the neutron and nuclear spin operators. 

The scattering probability can be separated into a 
factor depending on the initial and final neutron energy and 
a * scattering law^ , which depends only on the molecular 
dynamics of the scattering system. The expression on which 


U (r , R ^ * jR|^) 





the theory is based is given by 



Uer-e — — is the differential scattering cross-section 
dE^ 

per unit solid angle and per unit final neutron energy , for 


a neutron of tnomentum to scatter with momentum transfer 
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K = (2.3^) 

and energy transfer 

(k| - kj)/2m (2.35) 

S and S. are the scattering laws for coherent 
and incoherent scattering respectively. The above expression 
(equation (2.33) ) is for scattering from a system of N atoms 
which are identical except for a completely disordered dis- 
tribution of nuclear spins and isotopic species. The genera- 
lization to molecules with different atoms is straightforward. 
As mentioned earlier, t.e neutron spectra of materials contain- 
ing hydrogen, tx'.e most intense lines correspond to vibrations 
involving hydrogenous units and normally only the incoherent 
contribution from hydrogen is considered in analysing such 
spectra . 

For a general lattice, the differential scattering 
cross-section for incoherent scattering can be written as a 
series of terms giving the probability of absorption or 
emission of zero, one, or several phonons. 

<rinc 

d/LdE^ ^ dAdE^ 

where 1 is the number of phonons involved in the inter- 
action, The largest contribution to scattering is from 
the single-phonon (1=1 ) term with successive terms dimini- 
shing in magnitude. At low temperature, the main contribution 


dV 

^ xnc 


(2.36) 
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to the total incoherent scattering coines only from this 
"one-phonon" interaction and in practical calculations, 
higher order terras are neglected. This one-phonon appro- 
ximation contains the assumption that the vibrations of 
atoms about their equilibrium positions in the lattice are 
harmonic ana that phonon energies fico are greater than 


V' 


In the one -phonon approximation, the incoherent 

scattering cross-section for a polycrystalline sample can be 
written as; 


2 

(T, 


inc 


dA. dE, 


( 


^i PI exp (fiW/kgT - 1) 

exp g’^(o) 


Z-"i. 


xnc 


. . (2.37) 


where 0 is the scattering angle, is the incoherent 

th ±nc 

scattering length of the 1 nucleus in the unit cell, 

tii 

exp (-2¥^) is the Debye-Waller factor for the 1 nucleus 
and the sum over 1 includes all the n atoms in a unit 
cell. Here 


5 s(«- Wj(q)) ( 2 . 38 ) 

th 

where is the mass of the 1 nucleus, m is the 

neutron mass, V is the unit cell volume, q is the phonon 

wave-vector. The index j takes on 3n values corresponding 

•**1 — 

to 3ii branches of the dispersion relation. C.(q) are the 
polarization vectors mentioned previously. 
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In the case of monoatomic polycrystalline sample, the 
sum over 1 reduces to 


exp(_2¥^) exp(-2¥) g" ( UJ ) (2.39) 

where ^j^nc incoherent scattering length, m' is 

the mass of the atom and g ' ( tu ) is the effective frequency 
distribution, defined as in equation (2.38) above. The true 
phonon frequency distribution may be fornially defined as 





( 2 , 40 ) 


(2rr)^ ^ 

o 

Although g( to) and g" ( to ) are not directly 
related, comparison of equations (2.38) and (2,4o) shows 
that the true phonon frequency distribution has to be weighted 
by the square of the polarization vectors before it can be 
compared with the frequency distributions obtained from 
measurements of neutron scattering cross-sections. 

Although the frequency distributions obtained from 
equation (2.37) involve the use of several approximations, 
they do provide useful information and are found to be in 
qualitative, if not quantitative, agreement with calculated 
distributions. While relative intensities of the contribu- 
tions of the individual modes to g( to ) and shape of the 
individual peaks may be inexact, the frequencies of the 
modes can be obtained and identified. It should be noted 
that the cases where only the normal modes of vibrations 
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are obtained (not the dispersion relations) assuming the 
molecules to be free, equations (2,38) and (2*4o) can be 
approximately calculated without the integration over -the 
phonon wave-vector (q) . Relative intensity of the peaks 
in neutron spectrum helps in assignments of peak positions. 
The neutron spectrum of 1 ,3 }5-triamino-2 ,4 ,6-trinitrobenzene 
(TATNB) has been obtained and discussed in the next chapter. 

Infrared and Raman techniques used for the vibration 
spectra of TATNB (Chapter III) are very well developed^ ’ ^ ^ 
and are not discussed here. 
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CHAPTER - III 
VIBRATION SPECIRA OF 
1 , 3, 5 -IRIAMINO -2 ,4, 6 -TRINITROBENZENE 

1 , 3 1 5 -Triamino-2 ,4 , 6-trinitrobenzene (TATNB) is 
unusual in its thermal and solubility behaviour . It decom- 
poses above 300°C and dissolves best in concentrated sulphuric 
acid only. Its structure has been recently reported by Cady 
and Larson^ . The unit cell is triclinic with the space group 
pi and cell dimensions a=9»01A, b= 9*028 A, c = 6.812 A, 
and a = 1 08.59°, p = 91.82°, / = 119 . 97 °. There are two 
molecules per cell. It is found that the aromatic carbon- 

O 

carbon distance is 1 .44l A, the carbon-amino nitrogen distance 

, o 

is 1.316 A and there are six bifurcated nitrogen -oxygen 

hydrogen bonds. The hydrogen bonding is similar to the com- 

2 

pound 1 ,3-diamino -2 ,4 ,6-trinitrobenzene . The optical and 
cleavage properties of the crystal indicate that the molecules 
are arranged in planar sheets parallel to the a-b plane. This 
is confirmed by X-ray structure determination. The planar 
sheet structure greatly simplifies the calculation of normal 
vibrations of the system. The structure of the TATNB molecule 
is shown in Figure 2. Special interest in the molecule is 
because of the simultaneous substitution of electron donor 


and acceptor groups. 
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GROUP THEORETICAL ANALYSIS 

The molecule belongs to point group symmetry. 

There are 24 atoms in the molecule and 66 (24 x 3 - 6) funda- 
mental modes of vibration. The number of normal modes belong- 
ing to different symmetry species have been calculated using 
the relation (2.10). Character Table, optical activity and 
number of fundamental modes for each symmetry species are 
given in Table 1 . 

CALCULATIONS OF NORMAL VIBRATIONS 

A complete vibrational analysis has been carried out 
using Wilson GF matrix method^. The coordinates of the atoms 
are based on the average bond lengths: 

|(C-N)q = 1 .422 A, l(C-N)jj = 1.316 a, 

j/(C-C) = 1.441 A, ^(N-O) = 1 .242 A, 

|(N-il) = 0.925 A. 

All the angles between bonds are taken to be 120°. 

The internal coordinates used in thbse calculations are as 
f ollows : 

B ond s tr e t ching : 

r(C-N )q, r(N-O) , r(C-C), r(CHW)jj, r(N-H). 

In-plane angle bending : 

0(C-C-C)^^, 0(C-C~N)q, 0(C-N--O), 0(o-^r-o) , 0(C-C-C)jj, 

0(C— C— N)jj, 0(C— N—H) , 0(H— NHH) . 

Out -of -plane wagging ; 

6v>(C-N 02 ) , t*)(C-^H2 ) , W (R-N)q , O) )h • 



TABLE 1 


CHARACTER TABU), NUMBER OF NORMAL MODES FOR EACH 

SYMMETRY SPECIES AND OPTICAL ACTIVITY FOR 

— _______ 

SYMMETRY POINT GROUP 



E 

2C^ 

3C2 


2S3 

3 <rv 

a 

N 

Optical 

activity 


1 

1 

1 

1 

1 

1 

8 

Raman 

a ' 

^2 

1 

1 

-1 

1 

1 

-.1 

7 


e’ 

2 

-1 

0 

2 

-1 

0 

30 

Raman , I .R . 

A^ 

1 

1 

1 

_1 

-1 

-1 

2 


4 

1 

1 

-1 

-1 

-1 

1 

5 

I.R . 

E 

2 

^1 

0 

-2 

1 

0 

14 

Raman 


a 


’N* is the number of normal modes in the species 
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Torsional coordinates ; 

t(C-C), t(C-N)Q, t{C-^)jj 
The out-oT-plane modes include six benzene skeleton 
torsions and six torsions round (C— N) bonds. (C— C) torsion 

is sum of torsions due to (C— C) bonds and (C— N) bonds (sum 
of two torsions) whereas (C— N) torsion is sum of all four 
possible torsions, the resulting force constants should be 
accordingly interpreted. Using these definitions of internal 
coordinates G matrix is calculated and shown in Table 2. 
Symmetry coordinates have been written down as a normalized 
linear combination of internal coordinates. The definition 
of internal coordinates, and corresponding symmetry coordinates 
axe given in Appendix. Symmetry coordinates were used in 
construction of U-matrix which factors the G-matrix in a 
block diagonal structure. 

Urey-Bradley force field is used for the potential 

energy matrix. The expression for such a potential function 

4 

for the molecule is given by 


z 


K 


jh 




I- ■,,) 


K 


jk 


(A 




H'. „ r. . 
1 jk ij 


jk 


(A 


1 

2 


H 


i jk i j 


^ jk^^^ijk^ 


i,j,k 


^ ^ik ‘^ik '^ik^ 2 ^ik “^ik^ 

i, j,k 



V / A. n2 
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TABLJE 2 

G I4A.TRIX ELEMENTS FOR TATNB* 


I 

J 

G(I,J) 

I 

J 

G(I, J) 

1 

1 

0.155 

1 

2 

-0.036 

1 

24 

- 0.042 

1 

25 

0 . 1 CO 

1 

28 

-0.050 

1 

29 

0.100 

1 

33 

0.050 

1 

55 

-0.050 

2 

3 

-0.036 

2 

26 

-o.o 43 

2 

29 

-0.050 

2 

30 

-o.o 43 

3 

27 

- 0 .o 43 

3 

28 

-o.o 43 

4 

4 

0.167 

4 

5 

- 0.042 

4 

25 

-0.050 

4 

26 

-0.051 

4 

30 

-0.051 

4 

31 

-0.050 

4 

34 

- 0.055 

4 

36 

0.055 

4 

55 

0.050 

4 

56 

-0.050 

5 

7 

-0.036 

5 

8 

- 0 .o 42 

5 

31 

0.100 

5 

32 

-0.050 

5 

35 

0.134 

5 

36 

-0.067 

6 

6 

1 .063 

6 

7 

-0.036 

6 

34 

0 . 1 14 

6 

35 

-0.067 

7 

32 

o.o 47 

7 

33 

- 0 .o 47 

7 

36 

0.1 14 

8 

8 

0.167 

8 

25 

0.050 

8 

26 

-0.050 

8 

33 

0.105 

8 

34 

0.055 

8 

38 

0.101 

8 

39 

-0.051 

8 

43 

0.050 

8 

45 

-0.050 

9 

1 1 

-0.036 

9 

12 

- 0.042 

9 

37 

0.100 

9 

38 

-0.050 

9 

4 i 

0.100 

9 

42 

-0.050 

10 

10 

0.134 

10 

1 1 

-0.036 


I 

J 

g(i, J) 

I 

J 

G(I, J) 

1 

3 

-0.036 

1 

4 

- 0.042 

1 

26 

-0.050 

1 27 

-0.050 

1 

30 

-0.050 

1 31 

-O.05O 

1 

56 

0.050 

2 

2 

0.13^ 

2 

27 

0 .o 43 

2 28 

0.093 

3 

3 

0.134 

3 26 

o.o 43 

3 

29 

-0.050 

3 30 

0.093 

4 

8 

- 0.042 

4 24 

- 0 .o 42 

4 

27 

0.101 

4 28 

0.051 

4 

32 

0.105 

4 33 

- 0.055 

4 

37 

0.050 

4 39 

-0.050 

5 

5 

0.155 

5 

6 

-0.036 

5 

25 

-0.050 

5 26 

0.050 

5 

33 

-0.050 

5 34 

-0.067 

5 

37 

-0.050 

5 39 

0.050 

6 

32 

-o.o 47 

6 33 

o.o 47 

6 

36 

- 0 .o 47 

7 

7 

1 .063 

7 

34 

- 0.047 

7 35 

-0.067 

8 

9 

- 0 .o 42 

8 12 

- 0 .o 42 

8 

31 

- 0.050 

8 32 

- 0.055 

8 

36 

-0.055 

8 37 

-0.050 

8 

4 o 

-0.051 

8 42 

0.051 

9 

9 

0.155 

9 10 

-0.036 

9 

31 

-0.050 

9 32 

0.050 

9 

39 

-0.050 

9 ^0 

-0.050 

9 

43 

-0.050 

9 ^5 

0.050 

10 

38 

- 0 .o 43 

10 39 

0 .o 43 


. . Contd . 


*See footnote at the end of the table. 
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Table 2 (contd.) 


I 

J 

G(I,J) 

I 

J 

G{I 

:,j) 

10 

4 o 

0.093 

10 

4 l 

-0. 

050 

1 1 

38 

0.043 

1 1 

39 

-0. 

o43 

1 1 

42 

0.093 

12 

12 

0. 

167 

12 

31 

0.050 

12 

32 

- 0 . 

050 

12 

39 

0.101 

12 

4 o 

0 . 

051 

12 

44 

0.105 

12 

45 

- 0 . 

055 

12 

49 

0.050 

12 

51 

- 0 . 

050 

13 

15 

-0.036 

13 

16 

- 0 . 

042 

13 

43 

0.100 

13 

44 

- 0 . 

050 

13 

47 

0.134 

13 

48 

- 0 . 

067 

14 

14 

1 .063 

14 

15 

- 0 . 

.036 

14 

46 

0.114 

14 

47 

- 0 , 

.067 

15 

44 

0.o47 

15 

45 

- 0 , 

.047 

15 

48 

0.114 

16 

16 

0 , 

,167 

16 

37 

0.050 

16 

38 

- 0 . 

,050 

16 

45 

0.105 

16 

46 

0 . 

- 055 

16 

50 

0.101 

16 

51 

- 0 . 

.051 

16 

55 

0.050 

16 

57 

-0 

.050 

17 

19 

-0-036 

17 

20 

-0 

.042 

17 

49 

0.100 

17 

50 

-0 

.050 

17 

53 

0 

0 

0 

17 

54 

-0 

.050 

18 

18 

0.134 

18 

19 

-0 

.036 

18 

52 

0.093 

18 

53 

-0 

.050 

19 

50 

o.o43 

19 

51 

-0 

.043 

19 

54 

0.093 

20 

20 

0 

.167 

20 

25 

0.050 

20 

27 

-0 

.050 

20 

49 

-0.050 

20 

50 

-0 

.051 

20 

54 

-0.051 

20 

55 

-0 

.050 

20 

58 

- 0.055 

20 

60 

0 

.055 


I 

J 

G(I,J) 

I 

J 

0 ( 1 , J) 

10 

42 

- 0 .o 43 

1 1 

1 1 

0.134 

1 1 

4o 

- 0 .o 43 

1 1 

4 i 

-0.050 

12 

13 

- 0 . 042 

12 

16 

- 0.042 

12 

37 

-0.050 

12 

38 

-0.051 

12 

42 

-0.051 

12 

43 

-0.050 

12 

46 

- 0.055 

12 

48 

0.055 

13 

13 

0.155 

13 

l 4 

-0.036 

13 

37 

-0.050 

13 

38 

0-050 

13 

45 

-0.050 

13 

46 

-0.067 

13 

49 

-0.050 

13 

51 

0.050 

14 

44 

- 0 .o 47 

14 

45 

o.o 47 

14 

48 

- 0 .o 47 

15 

15 

1 .063 

15 

46 

- 0 .o 47 

15 

47 ■ 

-0.067 

16 

17 

- 0 .o 42 

16 

20 

- 0 .o 42 

16 43 

- 0.050 

16 

44 

- 0.055 

16 

48 

- 0.055 

16 

49 

-0.050 

16 

52 

-0.051 

16 

54 

0.051 

17 

17 

0.155 

17 

18 

-0.036 

17 

43 

-0.050 

17 

44 

0.050 

17 

51 

-0.050 

17 

52 

- 0.050 

17 

55 

-0.050 

17 

57 

0.050 

18 

50 

- 0 .o 43 

18 

51 

0 .o 43 

18 

54 

- 0 .o 43 

19 

19 

0.134 

19 

52 

- 0 .o 43 

19 

53 

-0.050 

20 

21 

- 0 .o 42 

20 

24 

- 0 .o 42 

20 43 

0.050 

20 

44 

-0.050 

20 

51 

0,101 

20 

52 

0.051 

20 56 

0.105 

20 

57 

- 0-055 

21 

21 

0.155 

21 

22 

-0.036 






Contd . 
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Table 2 (contd.) 


I 

J 

0(1, J) 

I 

J 

0(1, J) 

I 

J 

0(1, J) 

I 

J 

0(1, J) 

21 

23 

-0.036 

21 

24 

-0.042 

21 

25 

-0.050 

21 

27 

0.050 

21 

49 

-0.050 

21 

50 

0.050 

21 

55 

0.100 

21 

56 

-0.050 

21 

57 

-0.050 

21 

58 

-0.067 

21 

59 

0. 1 34 

21 

60 

-0.067 

22 

22 

1 .063 

22 

23 

-0.036 

22 

56 

-0.047 

22 

57 

0.o47 

22 

58 

0.1 14 

22 

59 

-0.067 

22 

60 

-0.047 

23 

23 

1 .063 

23 

56 

o.o47 

23 

57 

-0.o47 

23 

58 

-0.047 

23 

59 

-0.067 

23 

6o 

0.114 

24 

24 

0.167 

24 

25 

-0.050 

24 

26 

0,101 

24 

27 

-0.051 

24 

28 

-0.051 

24 

3^ 

0.051 

24 

31 

0.050 

24 

33 

-0.050 

24 

49 

0.050 

24 

50 

-0.050 

24 

55 

-0.050 

24 

56 

-0-055 

24 

57 

0.105 

24 

58 

0.055 

24 

60 

-0.055 

25 

25 

0.201 

25 

26 

-0.100 

25 

27 

-0.100 

25 

31 

-0.120 

25 

32 

-0.002 

25 

33 

0.122 

25 

34 

0.022 

25 

36 

-0.022 

25 

37 

0.020 

25 

39 

-0.020 

25 

49 

0.020 

25 

50 

-0.020 

25 

55 

-0. 120 

25 

56 

0.122 

25 

57 

-0.002 

25 

58 

-0.022 

25 

6o 

0.022 

26 

26 

0.197 

26 

27 

-0.097 

26 

28 

-0.117 

26 

30 

0.117 

26 

31 

0. 121 

26 

32 

0.002 

26 

33 

-0.123 

26 

34 

-0.022 

26 

36 

0.022 

26 

37 

-0.020 

26 

39 

0.020 

27 

27 

o^ 

• 

0 

27 

28 

0. 117 

27 

30 

-0.117 

27 

49 

-0.020 

27 

50 

0.020 

27 

55 

0.121 

27 

56 

-0-123 

27 

57 

0.002 

27 

58 

0 .022 

27 

60 

-0.022 

28 

28 

0.204 

28 

29 

-0.110 

28 

30 

-0.094 

28 

31 

-0.020 

28 

33 

0.020 

28 

55 

0.020 

28 

56 

-0.020 

29 

29 

0.220 

29 

30 

-0.110 

30 

30 

0.204 

30 

31 

0.020 

30 

33 

-0.020 

30 

55 

-0.020 

30 

56 

0,020 

31 

31 

0.201 

31 

32 

-0.100 

31 

33 

-0.100 

31 

37 

-0.120 

31 

39 

0.121 

31 

4o 

0.020 

31 

42 

-0.020 

31 

43 

0.020 

31 

45 

-0-020 

31 

55 

0.020 

31 

56 

-0.020 

32 

32 

0.213 

32 

33 

-0.113 

32 

34 

-0. l4l 

32 

36 

0.i4l 

32 

37 

0.122 

32 

39 

-0.123 

32 

4o 

-0.020 

32 

42 

0.020 

32 

43 

-0.020 


. . Contd. 
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Table 2 (contd.) 


I 

J 

0(1, J) 

I 

J 

G(I,J) 

32 

45 

0.020 

33 

33 

0.213 

33 

37 

-0.002 

33 

39 

0.002 

34 

34 

1 .391 

34 

35 

- 1 .285 

34 

39 

0.022 

35 

35 

2.569 

36 

37 

0.022 

36 

39 

-0.022 

37 

39 

-0.100 

37 

43 

-0.120 

37 

46 

0.022 

37 

48 

-0.022 

38 

38 

0.197 

38 

39 

- 0.097 

38 

43 

0.121 

38 

44 

0.002 

38 

48 

0.022 

38 

49 

-0.020 

39 

4 o 

0.117 

39 

42 

-0.117 

4 o 

4 l 

-0.1 10 

4 o 

42 

-0.094 

4 i 

4 i 

0.220 

41 

42 

-0.1 10 

42 

45 

-0.020 

43 

43 

0.201 

43 

49 

-0.120 

43 

50 

-0.000 

43 

54 

-0.020 

43 

55 

0.020 

44 

45 

-0.1 13 

44 

46 

- 0.141 

44 

51 

-0.123 

44 

52 

-0.020 

44 

57 

0.020 

45 

45 

0.213 

45 

49 

-0.002 

45 

51 

0.002 

46 

48 

-0.106 

46 

49 

-0.022 

47 

48 

-1 .285 

48 

48 

1 .391 

49 

49 

0.201 

49 

50 

-0.100 

49 

56 

-0.002 

49 

57 

0.122 

50 

50 

0.197 

50 

51 

- 0.097 

50 

55 

0.121 

50 

56 

0.002 

50 

60 

0.022 

51 

51 

0.197 

52 

52 

0.204 

52 

53 

-0.110 

52 

57 

0.020 

53 

53 

0.220 


I 

J 

G(I,J) 

I J 

G(I,J) 

33 

34 

0- i4i 

33 36 

-0. l4l 

33 

55 

-0.020 

33 56 

0-020 

34 

36 

-0.106 

34 37 

-0.022 

35 

36 

-1 .285 

36 36 

1 -391 

37 

37 

0.201 

37 38 

-0.100 

37 

44 

-0.002 

37 45 

0.122 

37 

49 

0.020 

37 51 

-0.020 

38 

4o 

-0.117 

38 42 

0.117 

38 

45 

-0.123 

38 46 

-0.022 

38 

51 

0.020 

39 39 

0.197 

39 43 

-0.000 

4o 4o 

0.204 

4o 

43 

-0-020 

4o 45 

0.020 

42 

42 

0.204 

42 43 

0.020 

43 

44 

-0.100 

43 45 

-0.100 

43 

51 

0.121 

43 52 

0.020 

43 

57 

-0.020 

44 44 

0.213 

44 

48 

0 . i4i 

44 49 

0.122 

44 

54 

0.020 

44 55 

-0.020 

45 

46 

0.i4i 

45 48 

-0.i4l 

46 

46 

1 .391 

46 47 

-1 .285 

46 

51 

0.022 

47 47 

2.569 

48 

49 

0.022 

48 51 

-0.022 

49 

51 

-0.100 

49 55 

0 

Cvi 

• 

0 

< 

49 

58 

0.022 

49 60 

-0.022 

50 

52 

-0.117 

50 54 

0.117 

50 

57 

-0.123 

50 58 

-0.022 

51 

52 

0.117 

51 54 

-0.117 

52 

54 

-0 - 094 

52 55 

-0-020 

53 

54 

-0.110 

54 54 

0-2o4 





. Contd . 
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Table 2 (contdo) 


I 

J 

0(1, J) 

I 

J 

G(I,J) 

54 

55 

0.020 

54 

57 

-0.020 

55 

57 

-0.100 

56 

56 

0.213 

56 

60 

0. l4l 

57 

57 

0.213 

58 

58 

1 .391 

58 

59 

-1 .285 

59 

60 

-1 .285 

60 

60 

1 .391 

6i 

68 

-0.04l 

61 

72 

-o.o4i 

6i 

77 

0.04? 

61 

78 

0.135 

62 

62 

2.977 

62 

67 

-0 . o44 

62 

73 

0.188 

62 

74 

-0.188 

62 

79 

-0.101 

62 

81 

0.101 

63 

69 

0.239 

63 

70 

-0.o4i 

63 

75 

-0.135 

63 

76 

-0.o47 

64 

64 

2.977 

64 

69 

-0-o44 

64 

74 

0.051 

64 

75 

0.188 

64 

81 

-0.101 

64 

83 

0.101 

65 

71 

0.239 

65 

72 

-o.o4i 

65 

77 

-0.135 

65 

78 

-0.o47 

66 

66 

2.977 

66 

67 

-0.o44 

66 

73 

-0.051 

66 

76 

0.051 

66 

79 

0.101 

66 

83 

-0.101 

67 

69 

o.o4o 

67 

71 

o.o4o 

67 

74 

-0.135 

67 

75 

0.046 

67 

78 

0.001 

67 

80 

0.279 

67 

84 

-0.279 

68 

68 

0.506 

68 

72 

o.o4o 

68 

73 

0-015 

68 

76 

0.046 

68 

77 

-0.046 

68 

81 

0.287 

68 

82 

-0.093 

69 

70 

-0.245 

69 

71 

o,o4o 

69 

75 

-0.001 

69 76 

-0.135 


I 

J 

G(I,J) 

I 

J 

G(I,J) 

55 

55 

0.201 

55 

56 

-O.IQO 

56 

57 

-0.113 

56 

58 

-0. l4l 

57 

58 

0.l4i 

57 

60 

-0. l4l 

58 

60 

-0.106 

59 

59 

2.569 

61 

61 

0.5o4- 

61 

67 

0,239 

61 

73 

-0-135 

61 

74 

-0.o47 

61 

80 

0.094 

61 

84 

-0 . 094 

62 

68 

0.294 

62 

69 

-0.044 

62 

75 

-0.051 

62 

78 

0.051 

63 

63 

0.5o4 

63 

68 

-0.04l 

63 

73 

0.047 

63 

74 

0.135 

63 

80 

-0.094 

63 

82 

0.094 

64 

70 

0.294 

64 

71 

-0.044 

64 

76 

00 

00 

• 

0 

1 

64 

77 

-0.051 

65 

65 

0.504 

65 

70 

-0.04l 

65 

75 

o.o47 

65 

76 

0.135 

65 

82 

-0.094 

65 

84 

0.094 

66 

71 

-o.o44 

66 

72 

0.294 

66 

77 

0-188 

66 

78 

-0.188 

67 

67 

0.480 

67 

68 

-0.245 

67 

72 

-0.24-5 

67 

73 

-0.001 

67 

76 

-0.046 

67 

77 

0.135 

67 

81 

-0.093 

67 

83 

0.093 

68 

69 

-0.245 

68 

70 

o.o4o 

68 

74 

-0.015 

68 

75 

-0.143 

68 

78 

0.143 

68 

79 

-0.287 

68 

84 

0.093 

69 

69 

0.480 

69 

73 

0.135 

69 

74 

0.001 

69 

77 

o.o46 

69 

78 

-0.046 





• m 

Contd . 
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Table 2 'contd.) 


I 

J 

0(1, J) 

I 

J 

Gfl , J) 

I 

J 

0(1, J) 

I 

J 

G(I , J) 

69 79 

0.093 

69 

80 

-0.279 

69 

82 

0.279 

69 

83 

-0.093 

70 

70 

0.506 

70 

71 

-0.245 

70 

72 

0.04o 

70 

73 

-o.o46 

70 7^ 

0-143 

70 

75 

0.015 

70 

76 

-0.015 

70 

77 

-0.143 

70 

78 

o.o46 

70 

80 

0.093 

70 

81 

-0.287 

70 

83 

0 .287 

70 

84 

-0.093 

71 

71 

0.480 

71 

72 

-0.245 

71 

73 

0.046 

71 

74 

-0 .o46 

71 

75 

0.135 

71 

76 

0.001 

71 

77 

-0.001 

71 

78 

-0.135 

71 

79 

-0.093 

71 

81 

0.093 

71 

82 

-0.279 

71 

84 

0.279 

72 

72 

0.506 

72 

73 

-0.143 

72 

74 

0.046 

72 

75 

-0.o46 

72 

76 

0.143 

72 

77 

0.015 

72 

78 

-0.015 

72 

79 

0.287 

72 

80 

-0.093 

72 

82 

0.093 

72 

83 

-0.287 

73 

73 

0 .2 10 

73 

74 

-0.054 

73 

75 

-0.006 

73 

76 

-0.107 

73 

77 

-0.006 

73 

78 

-0.037 

73 

79 

-0.117 

73 

80 

-0.108 

73 

81 

0.010 

73 

82 

0.107 

73 

83 

0. 107 

73 

84 

0.001 

74 

74 

0.210 

74 

75 

-0.037 

74 

76 

-0.006 

74 

77 

-0.107 

74 

78 

-0.006 

74 

79 

0.010 

74 

80 

-0.108 

74 

81 

-0.117 

74 

82 

0.001 

74 

83 

0.107 

74 

84 

0.107 

75 

75 

0.210 

75 

76 

-0.054 

75 

77 

-0.006 

75 

78 

-0.107 

75 

79 

0.107 

75 

80 

0.001 

75 

81 

-0.117 

75 

82 

-0,108 

75 

83 

0.010 

75 

84 

0.107 

76 

76 

0.210 

76 

77 

-0.037 

76 

78 

-0.006 

76 79 

0.107 

76 

80 

0. 107 

76 

81 

0.010 

76 

82 

-0.108 

76 

83 

-0.117 

76 

84 

0.001 

77 

77 

0.210 

77 

78 

-0.054 

77 

79 

0.010 

77 

80 

0.107 

77 

81 

0.107 

77 

82 

0.001 

77 

83 

-0.117 

77 

84 

-0.108 

78 

78 

0.210 

78 

79 

-0.117 

78 

80 

0.001 

78 

81 

0.107 

78 

82 

0.107 

78 

83 

0-010 

78 

84 

-0.108 

79 

79 

0.860 

79 

81 

-0.214 

79 

83 

-0.214 

80 

80 

12 .795 

80 

82 

-0.214 

80 

84 

-0.214 

81 

81 

0-860 

81 

83 

-0.214 

82 

82 

12.795 

82 

8^ 

-0.214 

83 

83 

0.860 

84 

84 

12 .795 











*Only nonzero elements are listed. Th.e numbering of th.e 
internal coordinates is in accordance with the order in 
which they are listed in the text. G matrix is symmetric, 
only elements on the main diagonal and above it are listed 
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•where /\r , /\0 . , At /\ Oj . are the internal coordi- 

nates corresponding to bond stretch, angle-bend, torsion and 
out-of-plane deformation respectively. The subscripts on the 
first two internal coordinates label the atoms involved. A^ t ^ 

■♦•Vi 

denotes the torsion about j bond and denotes the 

J 

j bond out-of -plane wag from the benzene-ring plane, 

A, q . , are the changes in the nonbounded distances (one-three 

4— -Ik. 

atom distances are taken for the purpose) . The primed quanti- 
ties are introduced as internal tensions and are taken as 
1/10 of the valence force constants. The Z - matrix defined 
by equation (2,5) was first constructed for TATNB . The F- 
matrix is calculated from Z -matrix and a set of initial 
force constants. These force constants were assumed, based 
on simpler molecules and they were then modified to get an 
agreement of the calculated frequencies with observed infrared, 
Raman and neutron frequencies. Force constants which give 
best fit are given in Table 3 and the corresponding F-matrix 
is shown in Table 4 . 

The problem of transferability of force constants has 

5 6 

been discussed in great detail by Scherer * . It is observed 

that the Urey-3radley force field can be successfully trans- 
ferred in some carbonyl, halogen containing molecules and 

7 

in s- trihalogen substituted benzenes • Similar successes 

were met with transfer of force field from symme-trical 

8 

perhaloethylenes to mixed chlorof luoroethylenes . 
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TABLE 3 

FORCE CONSTANTS (mdyn./A^ FOR TATNB^ 




K^(N - 

0 ) 

= 0.33 



k‘^(n - 

H) 

= 0.12 

K(C-N)^ 

= 1.98 



= 0. 12 

k(n-o) 

= 6.58 

K'^(R - 

n)h 

= 0.09 

k(c-c) 

= 3*86 




K(C-N)jj 

= 2.69 




k(nhh) 

= 5*87 

K*(C - 

c) 

= 0.22 



K*(C - 

">0 

= O-O 3 O 



K*(C - 

N)h 

= 0.0023 


H(C-C-C)q 

= 0,28 

F(C*C-C)^ 


0.30 

H(C-C-N)q 

= 0,88 

f(c*c*n)^ 


0.17 

H(C-i'J-O) 

= 1 .08 

F(C -N'O) 


0.23 

H(O-N-O) 

= 0.90 

F (0 ‘N'O) 


0.41 

H(C-C-C)p^ 

= 0.32 

f(c.c*c)j^ 

= 

0.31 

H(C-C-N)jj 

= 1 .14 

f(c-c*n)jj 

= 

0.39 

H(C-N-H) 

= 0.45 

F(C -N'H) 

= 

0.09 

H(H-N-H) 

= 0.35 

f(h*n*h) 


0.022 


stands for bond stretches, H for angle bends, for 

out-of -plane wags, for torsions and F for nonbounded 
interactions . 

V\ 

represents benzene ring. 
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TABLE ^ 

F- MATE IX ELEMENTS FOR TATNB* 


I 

J 

F(I, J) 

I 

J 

F(I,J) 

I 

J 

F(I,J) 

I 

J 

F(I,J) 

1 

1 

2.575 

1 

2 

0.178 

1 

3 

0.178 

1 

4 

0.132 

1 

24 

0.132 

1 

26 

0.094 

1 

27 

0.094 

1 

28 

0.124 

1 

30 

0. 124 

2 

2 

7.035 

2 

3 

0.318 

2 

29 

0.199 

2 

30 

0.113 

3 

3 

7.035 

3 

28 

0.113 

3 

29 

0.199 

k 

4 

4 .718 

4 

5 

0.302 

4 

8 

0 . 24 o 

4 

24 

0.232 

k 

25 

0.168 

4 

26 

0.095 

4 

31 

0-174 

4 

33 

0.215 

5 

5 

3.382 

5 

6 

0.069 

5 

7 

0.069 

5 

8 

0.302 

5 

32 

0.203 

5 

33 

0.203 

5 

34 

0-042 

5 

36 

0.042 

6 

6 

5.941 

6 

7 

0.017 

6 

35 

0.008 

6 

36 

0.033 

7 

7 

5.941 

7 

34 

0.033 

7 

35 

0.008 

8 

8 

4 .718 

8 

9 

0.132 

8 

12 

0.232 

8 

31 

0.174 

8 

32 

0.215 

8 

37 

0.168 

8 

39 

0.095 

9 

9 

2-575 

9 

10 

0.178 

9 

1 1 

0 

• 

00 

9 

12 

0.132 

9 

38 

0.094 

9 

39 

0.094 

9 

4 o 

0-124 

9 

42 

0.124 

10 

10 

7 .035 

10 

1 1 

0.318 

10 

4 i 

0.199 

10 

42 

0.113 

1 1 

1 1 

7.035 

1 1 

4 o 

0.113 

1 1 

4 l 

0-1 99 

12 

12 

4 .718 

12 

13 

0.302 

12 

16 

0.240 

12 

37 

0.168 

12 

38 

0.095 

12 

43 

0.174 

12 

45 

0.215 

13 

13 

3.382 

13 

14 

0.069 

13 

15 

0.069 

13 

16 

0.302 

13 

44 

0.203 

13 

45 

0.203 

13 

46 

0.042 

13 

48 

0.042 

i 4 

14 

5.941 

i 4 

15 

0.017 

14 

47 

0.008 

i 4 

48 

0.033 

15 

15 

5.941 

15 

46 

0.033 

15 

47 

0.008 

16 

16 

4.718 

l 6 

17 

0.132 

16 

20 

0.232 

16 

43 

0.174 

16 

44 

0.215 

1 6 49 

0.168 

16 

51 

0.095 

17 

17 

2.575 

17 

18 

0.178 

17 

19 

0.178 

17 

20 

0.132 

17 

50 

0.094 

17 

51 

0.09^ 

17 

52 

0.124 

17 

54 

0.124 

18 

18 

7.035 

18 

19 

0.318 

18 

53 

0.199 

18 

54 

0.113 

19 

19 

7.035 

19 

52 

0-1 13 

19 

53 

0.199 

20 

20 

4.718 

20 

21 

0.302 

20 

24 

0.240 

20 49 

0.168 

20 

50 

0-095 

20 

55 

0.174 

20 

57 

0.215 


. . Contd. 

*See footnote at the end of the table. 



Table 4 (contd.) 


I J F(I,J) I J F(I,J) I J F(I,J) I J F(I,J) 


21 21 
21 56 
22 22 

23 23 

24 25 

25 25 
29 29 

33 33 
37 37 
4 i 4 i 

45 45 
49 49 
53 53 
57 57 

6i 6l 

65 65 
69 69 
73 73 
77 77 
81 81 


3.382 

0.203 

5.941 

5.9^1 

0.168 
0.482 

1.106 
1 .380 

0.482 

1.106 
1 .380 
0.482 

1.106 
1 .380 
0.330 
0.330 
0.119 
0.225 
0.225 
0.030 


21 

22 

21 

57 

22 

23 

23 

58 

24 

27 

26 

26 

30 

30 

34 

34 

38 

38 

42 

42 

46 

46 

50 

50 

54 

5 ^ 

58 

58 

62 

62 

66 

66 

70 

70 

74 

74 

78 

78 

82 

82 


0.069 

0.203 
0.017 
0.033 
0.095 
0.993 
1 .213 
0.480 

0.993 
1 .213 
0.480 

0.993 
1 .213 
0.480 
0. 124 
0.124 
0.090 
0.225 
0.225 
0.002 


21 

23 

21 

58 

22 

59 

23 

59 

24 

55 

27 

27 

31 

31 

35 

35 

39 

39 

43 

43 

47 

47 

51 

51 

55 

55 

59 

59 

63 

63 

67 

67 

71 

71 

75 

75 

79 

79 

83 

83 


0.069 

0.042 

0.008 

0.008 

0.17^ 

0.993 

0.529 

0.356 

0.993 

0.529 

0.356 

0.993 

0.529 

0.356 

0.330 

0.119 

0.119 

0.225 

0.030 

0.030 


21 

24 

21 

60 

22 

60 

24 

24 

24 

56 

28 

28 

32 

32 

36 

36 

4 o 

4 o 

44 

44 

48 

48 

52 

52 

56 

56 

60 

60 

64 

64 

68 

68 

72 

72 

76 

76 

80 

80 

84 

84 


0.302 

0.042 

0.033 
4.718 
0.215 
1 .213 
1 .380 
0.480 
1 .213 
1 ,380 
0.480 
1 .213 
1 .380 
0.480 
0.124 

0.090 

0.090 

0.225 

0.002 

0.002 


*0n.ly nonzero elements are listed. The numbering of the 
internal coordinates is in accordance with the order in 
which they are listed in the text (p. 29 ). F matrix is 
symmetric , only elements on the main diagonal and above 
it are listed . 
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EXPER IMENTAL 

The infraref' spectrum of polycrystalline sample of 
TATNB in the form of KBr pellets were obtained using Perkin 
Elmer 521 spectrometer in the region 5000 - 250 cm ^ and is 
shown in Figures 3 and 4. Laser-excited Raman spectrum has 
been recorded on Coderg spectrophotometer and shown in 
Figure 5- The inelastic neutron spectrum of TATNB recorded 
at the AMMRC slow chopper is shown in Figure 6. This was 
also helpful in assignments of the frequencies. Since TATNB 
is insoluble in almost all the solvents it was not possible 
to get an infrarec or Raman spectrum in solution form. 

This made deuteration impossible. 

DISCUSSION 

The calculated and observed frecpencies in I .R . , Raman 
and neutron scattering and their possible assignments based 
on potential energy distribution, symmetry principles, peak 
intensities, bond order and group frequencies are given in 
Table 5. The detailed discussion for assignments of the 
modes follows. 

NH^ Group Frequencies 

The intense broad bands in infrared spectrum at 3328 
and 3227 cm~^ are assigned to NH 2 asymmetric and symmetric 
stretches respectively. These modes appear at 3^^0 and 
33^0 cm~^ in aniline^ and at 3^70 and 3350 cm in 
p-nitroaniline^®. The shift towards lower frequencies in 
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Fig. 5 Raman spectrum of TATNB 
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Fig. 6 Inelastic neutron spectrum of TATNB , j 
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TABLE 5 



CALCULATED AND 

OBSERVED frequencies FOR TATNB^ 



{ 

a) In -plane 

modes 

Species 

Calcu- 

lated 

Obser- 

ved 

Experi- 

ment 

Assignment^ 

a; 

3227 


- 

NHg Sy. Stretch 


1618 

- 

- 

NH^ Sy. Bending 


1250 

1303 

Raman 

NOg Sy. Stretch 


1 l 4 l 

1154 

Raman 

X - Sens. (C-NH^) Stretch 


1017 

1017 

Raman 

X - Sens. (C— NO^) Stretch 


678 

687 

Raman 

NOg Sy. Bending 


476 

510 

Raman 

X - Sens. 0 (C-C— C)jj 


243 

236 

Raman 

X - Sens. 0 (C— C— C)q 

\ 

^2 

3329 

— 


NHg Asy. Stretch 


1544 

- 

- 

NOg Asy. Stretch 


l 4 l 1 

• 

- 

Ring Stretching 


1130 

1 129 

Raman 

NHg Rocking 


823 

- 

- 

0 (C— C— N)jj + 0 (C— N-H) 


421 

430 

Neutron 

NO^ Rocking 


203 



X - Sens. 0 (C— C— N)q 

. . (Contd .) 


(See footnote at the end of idie table). 
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Table 5 (contd.) 


(a) In -plane modes 


Species 

Calcu- 

lated 

Obser- 

ved 

Experi- 

ment 

Assignment^ 

1 

E 

3329 

3328 

I.R . 

NHg Asy . Stretch. 


3^27 

3327 

I,R . 

NHg Sy. Stretch 


I6I8 

1618 

I.R . 

Asy. Bending 


1572 

1571 

I.R . 

NOg Asy. Stretch 


1468 

1447 

I.R. 

Ring Stretching 


1312 

1322 

I.R . 

Ring Stretching 


1228 

1229 

I .R . , Raman 

NO2 Sy. Stretching 


1029 

1033 

I.R . 

NH2 Rocking 


907 

880 

Raman , I .R . 

X - Sens. (C-NHg) Stretch 


731 

725 

I.R . 

NO2 Bending 


589 

606 

Raman , I .R . 

X - Sens. (C--NO2) Stretdi 


447 

445 

I .R . , Raman 

NO2 Rocking 


354 

360 

Neutron, Raman X - Sens. 0 (C— C— 


24 o 

- 

- 

X - Sens. 0 (C-C-C)q 


126 

125 

? Raman 

X - Sens. 0 (C-C-N)q 


, , (Contd,) 


^9 


Table 5 (contd.) 

(b) Out -of -plane modes 


Species Calcu- 
lated 

Obser 

ved 

Experi- 

ment 

A s sxgnment^ 

tt 

219 

220 

Neutron 

(C— NH^) Torsion 


148 


- 

(C-NO^) Torsion 

11 

•^2 

799 

782 

I.R. 

NH^ Out of plane wagging 


572 

567 

I.R. 

NO 2 Out of plane wagging 


507 

510 

Neutron,! .R . 

Ring distorsion 


I89 

- 

- 

(C-^) Torsion 


74 

- 

- 

Ring distorsion 

11 

E 

805 

8I9 

Raman 

NH^ Out of plane wagging 


601 

606 

Raman 

NO^ Out of plane wagging 


388 

00 

Raman 

Ring distorsion 


302 

- 

- 

Ring distorsion 


215 

212 

Raman 

(C-NH^) Torsion 


169 

173 

Raman 

t(C-C) + t{C--N)Q 


1 17 

125 

? Raman 

t(C— NO^) + t(C— C) or 





Lattice Mode 

a 



-1 


All 

frequencies are 

xn cm ♦ 


^ A 

0 » 

deformation 

mode 5 

t, torsional 

mode ; Sy . , Symmetric ; 


Aay . , Asymmetric ; Sens . , 


Sensitive . 
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TATNB is due to (a) intra and inter molecular hydrog^en 

bonding and (b) due to resonance character of the amino 

group. Resonance forms of TATNB are shown in Figure 7- 

The empirical relations between the V* (NH„) and i;^(NH^) 

a. s ^ s ^ 

1 1 

modes for amines do not hold for TATNB. The obvious 
reason is due to large perturbations and crowding of mole- 
cules in the crystal. However, the broadening of half- 

12 

widths is in conformity with hydrogen bonding . In case 
of bending and other modes the shift due to hydrogen bonding 
is in the opposite direction. The H— N— H angle bend appears 
at l6l8 cm~^ and the rocking, wagging and torsional modes 
are at 1033* 782 and 220 cm”^ respectively. The last two 
assignments are justified from the neutron spectrum in which 
the intensity of a mode is weighted according to the amplitude 
of motion of the hydrogen atoms. One would expect larger 
displacements in wagging rather than torsional modes. Hence 
the peak at 760 cm“^ in the neutron spectrum is more likely 
due to NHg wagging motion and one at 220 cm due to torsional 
motion. The peak at 760 cm ^ appears broad due to the 
anharmonicity arising from hydrogen bonding of the protons. 
This mode is also seen in I.R. spectrum at 782 cm . In 
aniline^ NHg wagging, rocking and torsional modes appear at 
670, 1050 and 220 cm"^ respectively. NHg rocking (A^ mode) 

* y : Asjwmetric stretch . 5 Synanietric stretch# 






5Z 


appears in calculations at 1 1 30 cm“^ and the weak Raman band 
at 1129 cm ’’ should correspond to this mode. Though, 
species is Raman inactive, the appearance of this mode 
could be due to departure from complete symmetry of the 
molecule . 

NOg Group Frequencies 

NO^ group symmetric and antisymmetric stretch frequen- 
cies have been observed in the narrow range of ( 1338 - 1355 ) cm- 
and (151O-I555) cm respectively in para substituted nitro- 
benzene derivatives^^. The Raman band at 1302 cm~^ should 
correspond to totally symmetric NO2 stretch mode. All the 
totally symmetric modes are expected to appear strong in Raman 
spectrum. A shift towards lower frequencies is a result of 
strong hydrogen bonding in the molecule. Such a situation 
is also seen in case of hydrogen bonded 0 -nitrophenol where, 
falls to 1320 cm , A study of a series of multiply- 
substituted nitro -benzene compounds shows that the bands 
corresponding to (NOg) stretch appear in the I .R , spectrum 
with fairly large intensity. Two intense bands at 1229 and 
1181 cm ^ may be because of degenerate (NO2) mode- 

Splitting of the degenerate mode arises because the effective 
symmetry of a single constituent molecule in the crystal may 
be lowered on account of its environment and selection rules 
for free molecule may no longer be obeyed. Departure from 
the normal range of the asymmetric (NO2) stretch mode is 
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either as a result of heavy multiple substitution or steric 

effects. In 1 , 3 > 5-trinitr obenzene this mode appears at 
_ 1 

1557 cm while 1 ,2 , 3-trinitrobenzene has two bands at 1558 
_1 13 

and 1572 cm . The non-degenerate and degenerate 

-1 

bands in our calculations are obtained at 1544 and 1572 cm 

respectively. However, only the latter mode is observed in 

I .R . and has a large mixing in potential energy with (C— C) 

stretch (as seen from the calculations). Hydrogen bonding 

has little effect on < (N0„) band and even in 0-nitrophenol 

as ' 2 ' 

the U band is at 1537 cm""^ . 

3. S 

There is a controversy about the assignment of (NO 2 ) 

14 

symmetric bending of nitrobenzene. Green et al , have 

-1 15 

assigned this mode at 677 cm” while Stephenson et al. 

_ 1 

assignment corresponds to a mode at 850 cm . In 1 ,3j5- 

*1 ^ 

trinitrobenzene the mode appears in Raman spectrum at 
714 cm”^ and favours the assignment of Green at al. The band 
at 687 cm”^ in Raman scattering of TATNB has been assigned 
to totally symmetric (NO 2 ) bending mode. The absorption 
bands in I .R . spectrum at 701 and 725 cm”^ corresponds to 
the degenerate (NO 2 ) bending mode. The shift towards somewhat 
higher frequencies is due to hydrogen bonding. The out-of- 
plane wagging mode of NO 2 group is expected to have higher 
frequency than in-plane mode because (a) the moment of 
inertia of the group about an axis perpendicular to its plane 
is much greater than that about an axis in its own plane 
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ana (b) the repulsion in F.E. function for the out -of -plane 
mode is probably greater than that of the in-plane mode. 

Hence the band at 445 cm~^ is assigned to the in-plane rocking 
mode (E* ) , this band appears in Raman shifts at 444 cm~\ A 
weak I.R. mode at 5^9 cm ^ is assigned to NOg out-of -plane 
wagging which is confirmed by the assignments of the nitro- 
benzene at 532 cm \ In 1 ,3 »5-trinitrobenzene^'^ the in-plane 
rocking and out-of -plane wagging modes appear at 450 and 
518 cm ^ respectively. The corresponding degenerate mode {^^(NOg) 
could be Raman band at 606 cm~ . ^^2 mode appears 

-1 15 

at 139 cm in nitrobenzene and oxor calculations give at 

148 cm ^ (a!| species) . Some alternative assignments have 

* 1 3 ^0 

been suggested for NOg torsional mode ~ . But still it 

seems that no definite assignment of this mode has been made. 
X-Sensitive Modes 

1 Stretching Modes 

A study of symmetric -trihalogen substituted benzene 
compounds^ shows that x-sensitive stretching modes lie in 
a broad range and wheii these compounds are deuterated 
completely they shift towards lower frequencies. In TATNB 
since all the hydrogens are replaced by NH^ and NOg group 
symmetrically, we expect the x-sensitive frequencies to lie 
in the lower frequency side of the above range. However, 
one should be cautious in comparing these frequencies with 
normal (C—N) stretches, since in such cases the stretching 
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ot (C— N) bond will necessarily involve ring stretching and 
motion of hydrogens or oxygens as the case may be. The 
sharp Raman band at 1154 cm"^ has been assigned to (C— NHg) 
stretching mode ( ) . All the totally symmetric modes are 
expected to appear strongly in Raman spectrum. A very large 
intensity in this band is because of large degree of covalent 
character. Since the change in polarizability during stretch- 
ing is related to the extent to which electrons are shared in 
the bond i.e. the bond order, (C— NO^) x- sensitive stretching 

mode will be expected to be lower than (C— NH^) mode and as 

- 1 

such has been assigned at 1019 cm . The potential energy 

contributions in these modes are only 48'^ and 19 ^ from 

(C-NHg) stretch and (C-NO^) stretch respectively. In the 

1 6 

Raman spectrum of 1 , 3 i 5 -tr±nitrobenz;ene the (C— NO 2 ) stretch 

mode has not been observed because of its low intensity* The 

- 1 

corr esponding degenerate modes could be 88 O and 606 cm from 
Raman shifts, supported by calculations. 

2 Bending and Deformation Modes 
The neutron peaks at ^30 and 360 have been assigned 

to (C--CHM) deformation modes on the basis of shift in the 
Raman scattering, (C— C— N)jj bending frequency in aniline and 
normal mode calculations. The (C-C-C)jj^q bending modes have 
been designated as x- sensitive since they involve considerable 
(C— N) stretching. They have been assigned at 510 cm and 
236 cm”^ from Raman scattering. The modes are in expected 
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region as seen from nitrobenzene^^, aniline'^ and p-nitro- 

aniline . The sharp though weak peak in the neutron spectrum 
- 1 

at 510 cm could arise due to either out-of-plane wagging 

motion of the nitro-groups or the ring deformation arising 

from the out -of -plane motions of the nitrogens of the nitro 
, - 14 

and amxno group , The latter is more likely to involve the 
displacement of protons and hence expected to appear in 
inelastic neutron scattering. 

Ring Modes 

The ring stretching modes (e’ species) have been 
assigned at 1447 and 13S9i cm ^ from I.R, spectrum. A consi- 
derable lowering of these modes from the benzene stretching 
modes is as expected since there is a large increase in (C— C) 
bond length in TATNB (l.44l 1 ). A mode at 101? cm“^ from 
Raman scattering corresponds to ring breathing mode of 
benzene, but as this mode involves considerable (C-N) stretch, 
it has been classified as x-sensitive stretching mode discussed 
above. The out-of -plane ring distortion mode is expected 
below 250 cm~^ which is the limit of Perkin Elmer equipment. 

The peak at 110 cm”^ and broad band with a maximum at 
65 cm~^ in neutron spectrum are in the region of lattice modes 
involving the stretching and bending of hydrogen bonds and 
spread nearly over the region of lattice modes observed in 
other hydrogen bonded systems. Raman shifts also support 


this - 
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Diagonal elements of the mean square amplitude matrix 
have been calculated using relation (2.23). A spatial repre- 
sentation of some of the normal modes are shown in Figures 8 
and 9 * The thermodynamic parameters i.e, entropy, heat content, 
heat capacity and free energy for TATNB molecule are obtained 
from relations ( 2 . 24 ) to ( 2 . 29 ). They are given in Table 6. 

At the end, it should be added that calculations on an 
isolated molecule is a sin^lified model when one considers 
actual conditions prevailing in solid state. This is especially 
true for TATNB where hydrogen bonds play an important part. 

The effects that would be expected if one considers the 
symmetry of the space group rather than point group, are dis- 
cussed in the VI Chapter. The work on isolated molecule should, 
however , be a good starting point . 
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TABLE 6 


THERMODYNAMIC QUANTITIES FOR TATNB^ 


Temp . 

eA) 

Heat 

Content 

( o’ 

\ T / 

Heat 

Capacity 

(CP°) 

Free 

Energy 

f F°-E° 1 

0 \ 

1 ! „ 

\ T ^ 

Entropy 

(sc) 

300.0 

39-62 

68,50 

88. 4o 

128 .02 

400.0 

48,51 

81 .22 

101 .05 

149 .56 

500.0 

56-10 

91 .10 

112.71 

168 .80 

573.0 

60.95 

96.85 

120.68 

181 .62 

6oo.o 

62.61 

98 -72 

123 .52 

186.13 

700.0 

68.22 

104.65 

133.61 

201 .82 

800.0 

73.09. 

109.37 

i 43 .o 4 

216.13 

900.0 

77.35 

1 13 .20 

151-90 

229.25 

1000.0 

81 .11 

116.36 

160.25 

24i .35 


®'A11 quantities are in cal./™ol-/° Absolute. 


9 
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CHAPTER - IV 

SMALL-ANGLE X-RAY SCATTERING - THEORETICAL DETAILS 

The folded chain lamellae with a rather uniform thick- 

O 

ness L of between one and few hundred A and much larger 
lateral extension are the main structural unit of the polyn^r 
solid. They form spherulites originating from primary crysta- 
llization nuclei which, by secondary nucleation and to a large 
extent by non-crystallographic branching, grow first into a 
bundle-like and later into a radial arrangement of ribbon- 
like, more or loss twisted, lamellae. The size and perfection 
of spherulitic structure depend on crystallization conditions, 
on temperature, density of nuclei, impurity concentration, and 
so on. The parallel stacking of lamellae creates a. nearly 
ideal periodic lattice which, due to large period L, gives 
rise to small-angle X-ray scattering. The intensity of 
scattering depends on the square of the electron density 
difference between the crystal lattice and the region between 
the adjacent crystals, which contains chain folds, tie mole- 
cules and free ends of chains incorporated only partially 
in the crystalline lattice. The thickness of the density- 
deficient layer is the second factor affecting the scattering 
intensity. Small-angle X-ray scattering, however, does not 
yield any information about the fine structure of the inter- 
crystal layers, particularly whether they are amorphous or 
and how long the folds are and so on. 


not , how many 
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For polyethylene, layer like structure was first 

indicated in 195^ by Claver et al . who compeired the free 

surface of a bulk sample to a fanned deck of cards. Later, 

the occurrence of layers was directly confirmed by electron- 

2 

micrographs, both of fracture surfaces and of ultramicrotome 

3 

sections and also by elec tr on- micro sc opic examination of the 

debris obtained by oxidation of polyethylene with fuming 

. 4 

nitric acid . A .detailed analysis of small-angle X-ray sca- 
ttering curve can reveal information about the arrangement 
of crystallino and amorphous regions within the layers - 

A reasonable-model calculations can be performed to 
explain the scattering curve. If the model is correct, the 
parameters of the model can be varied to make the calculated 
scattering curve fit the experimental one. However, for the 
studios of semicrystalline materials it appears advantageous 
to calculate the Fourier transform of the scattering curve 
as an intermediate step. The Fourier transform of the sca- 
ttering curve gives the position-correlation fimction as 
introduced by Debye and Bueche^. The correlation function is 
proportional to the autoconvolution of electron density var- 
iations in the sample and can be evaluated for an assumed 
model. The calculated correlation function can be compared 
with that derived from the experimental scattering exarve and 
the variables of the model can be best adjusted. 
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EXPERIMENTAL CORRELATION FUNCTION 

The intensity of small-angle scattering from a single 
region of parallel layers along normal direction is represen- 
ted by an intensity function i(s) in reciprocal space. This 
is given by 


i(s) = 2 V j y (x) Cos 2 fjxs dx 

J 


(4.1) 


where 


y (x) = J ^ ^ 


(4.2) 


and s = 


o 

2 Sin 0 


V = Irradiated volume of the region under considera- 
tion 

X = Coordinate perpendicular to the layer 
9 = Half-angle of diffraction 

^ = Wave-length of X-rays 

^( ^ ) = Local fluctuations of electron density around 
the average value. 

From equation (4.1), we get 


f(D 


/ (x) = ^ J i(s) Cos ZTTxs ds 


t ^ roo 

and y'(0) = <1^ > Ty J 


ds 


(^. 3 ) 

(4.4) 


is the average of the square of the electron 


density fluctuations. 


The correlation function as defined by Debye and Bueche' 


is f 
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n^) = 


■y. \ 

’ ' zy 




^(o)^ ’ ■'<>'(^)<+ 1) (4.5) 


When all orientations of the sets of layers are present, 
the measured intensity l(s) has spherical symmetry and 
is related to i(s) by 


I(s) i(s) / 4 TT (4,6) 

The relation between )^(x) eind the experimental 
scattering curve obtained from relations (4.3) to (4.6) is 
given by 


roo 2 

J S l(s) Cos (2 TTxs) ds 



o 


(4.7) 


Using the experimental intensities l(s), the position 
correlation function can be obtained from the relation (4.7). 
Some difficulties are involved in evaluation of this equation. 
Experimental scattering curve can only be obtained from a 
lower limit of the scattering angle. The full small-angle 
scattering curve is obtained by extrapolating to zero -angles 
the observed curve. It is foimd that the method of extra- 
polation has very little effect on the relevant part of the 
correlation function, neither the position nor the height 
of the first maximum is seriously affected, not even if 
absurd methods of extrapolation are used. However, the tail 
of the scattering curve does affect the correlation curve. 
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specially at small values of x. In order to correct for 
the overlap of the small and large angle scattering curves, 
the base line is assumed tangential to the tail end of the 
experimental curve. The part of the correlation curve 
affected by this procedure is shown dotted (V Chapter - 
Figures 17 and 18). 


THE CORRELATION FUNCTION FROM A MODEL 


A layer structure model with alternate high and low 
electron density is assumed. These layers are designated 
as crystalline and amorjhous, though the use of these names 
do not imply that the crystalline layers are entirely crysta- 
lline or the amorphous layers are entirely amorphous. The 
thicknesses of these layers are assumed to have a Gaussian 
distribution. If x and x^ be the thicknesses of crysta- 
llxno and amorphous regions then normalized distribution 
functions around mean values c and a are given by 


P (x ) 
c '■ c ' 


P (x ) 

a' 


1 ) 


(4.8) 




1 

IT 


(^) 


(^•9) 


y277 

jj, c and are the widths of these functions. 

^CO 

Average crystallite length, C = j 


. (4.10) 



j roo 

i X P fx > dx 

a a ^ a ^ ' 


Crystallinity (0) is defined as 


0 


C +A 


(4.12) 


Following’ the method of Debye, Anderson and Brumberger 
the position correlation function for a two-phase model is 
given by the expression 


,, Pc>>- 0 

y(^) = (^'•‘3) 

p (x) denotes the probability that if one end of a 

o c 

hypothetical measuring rod of thickness x is in the crysta- 
lline region (higher electron density phase) then the other 
end is also in one of the crystalline regions . The situation 

is very much simplified by the fact that the polyethylene 

7-9 

chains are found to be normal to the layers Let 

PccC^c) bo expressed as 

Pec %ac + %acac "" 

whore q (x) denotes the probability of finding the two 

ends of the measuring rod in the same crystalline region, 

q (x) to that of the ends being in neighbouring crysta- 
cac ^ 

nine regions and so on. 

The probability (x) for a length x in a 

crystalline region of thickness x^ is " ^)/^c 



seen from Figure 10(a). For 


ail the thicknesses we get 


%oo 


q (x) 


J - X) Pjx^) dx 


fOO 

^ X P (x ) dx 

C ri ' n ' r 


(4x15) 


'"c C'' C' ■ c 


I' or contributions to measuring rod must 

be Situated as shown in Figure lo(b). In this case, it is 
composed of throe vectors x^ , x^, x^ . According to a well 

J 

known theorem , the distribution of such Vasum of vectors 
(Vector -Sum) is given by the normalized distribution function 




(4.16) 


Sign N denotes the convolution product of the two 
functions. given by the equation (4.9) and 

contributions for Q(3 c^) are obtained 
from tho crystalline regions which has thicknesses more than 
x^ i.e. tho fraction: 

/ on 

P (x ) dx 
o'- c' c 

Normalizing this function, we get 



(4.17) 


Using relations (4,l6), (4.17), (4,8) to (4,11), the 
expression for P ^^(x) can be written as 

C8iC ^ 





Now, tho relation between q (x) and 

C cLC ^ 

will bo dovoioped. The chain length can take only N 


discrete values. q (x) can be written as 

xac ^ ^ 


q ( 3c) 

xac ^ ' 


^cac(^) 


( 4 . 13 ) 


T (x'i is tho number of positions of a measuring 
cac' ’ 

rod of length x with A and B in neighbouring crysta- 
lline layers and R (x) is the number of positions of a 

c 

measuring rod when A is in crystalline layer and B 
in either of layers. R^(x) is equal to 0N per unit 

length. S , , is the total number of vectors of type cac . 

C «iC 

This number is equal to 0^N^(C+A) per unit length , T^ac^*^ 
is given by •' 


T 

cac 


(*) 


^cac^^^ ^cac 
N 


( 4 . 20 ) 


putting the respective values we get 


q 

xac 


(x) = 0(C+A) 


(4.21) 
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%acac(*) = Pcacac 

Pcaoac(*) = CPa^r^a'^ ('*•^3) 

A relation similar to equation (4.18) can be written 
for P (x) and higher order terms. 

The expression for correlation function is obtained 
using relations (4,13) to (4,15)? (^*21) and (4.22) and is 
given below : 

a roo 

y(^'l = I 5 / ^oac 

^ i '^X 

* + - 0 j 

(4.24) 

The theoretical correlation curve equation (4.24) 
has been numerically solved. Theoretical curve has been 
shown in Table 7 (Chapter V) for some particular values of 
model parameters (c, a, /j-c, A a) . These parameters can 
bo varied to get an agreement between theoretical and 
experimental correlation curves. 


Similarly 

wher e 
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CHAPTER - V 

■§.?1^L-ANGLE X-RAY SCATTERING FROM POLYETHTLENE AM) WOOL 

Morphological structure of polyethylene is a subject 
of considerable interest because of its iirqportance in under- 
standing the process of crystallization and thermodynamic 
behaviour . It is now almost established that polyethylene 
has layer like structure^ but still there is much contro- 
versy about the detailed arrangement of the crystalline and 
amorphous regions in these layers. As seen in the last 
chapter , small-angle X-ray scattering depends on the differ- 
ence of the square of the average electron densities between 
different regions. The analysis of such a scattering will 
give information regarding the detailed arrangement of the 
layers. Abundant experimental observations have been reported 

earlier for melt-crystallized and fibrous polyethylene 
9-24 

samples , Howover , the experimental observations on bulk- 

crystallized samples are much less. Secondly, in most of 

the studies the interpretation of the scattering maxima has 

been done assuming the Bragg condition to be valid. This is 

obviously not the right approach, since normally the arrange- . 

ment of ordered and disordered regions is not regular. By 

23 

way of exanaple, the results obtained by Geil from small- 

angle X-ray scattering are not in agreement with their 

24 

electron-microscopic measurements. Recently, Vonk et al. 
have reported a study of correlation functions based on a 
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two-phase parallel-layer model with different average electron 
densities. However, their results vary within wide limits. 

25 

Russian workers have reported results in agreement with 
those obtained from electron-microscope. In view of these 
discrepancies, it was decided to study the position correla- 
tions in a linear polyethylene sample which has been annealed 
just below the melting point and cooled very slowly to room 
temperature. To test the validity of the two-phase model the 
experiments are also performed on irradiated polyethylene 
sample . 

EX PERIMENTAL 

26 

The Luzzati-Baro small-angle scattering camera has 

been used with linear collimation system. The three slit 

system is used to remove pairasitic scattering. The slit 

system is shown in Figure 1 1 . The camera is evacuated to 

avoid air scattering. The alignment corresponds to a Bragg 

value of about 300 A for CuK radiation. Photographic tech- 

ot 

nique is used for recording the scattering pattern- Sample 
to film distance was 10. 8 cms, for polyethylene and 11.4 cms. 
for wool satcples. CuK radiation with nickel filters was 
used. Intensity measurements were made on a Kipp-Zonan 
microphotometer. The intensities were converted to relative 
scale with the help of a log-scale obtained by means of a 
smooth wedge. Polyethylene samples were about 1 mm thick. 

The X-ray patterns of Polyethylene and Irradiated Polyethylene 




STEM 
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sairples and a log-scale variation are shown in Plate 1 
(a,b,d) and their microphotometer records are shown in 
Figures 12, 13 and respectively. 

SLIT CORRECTION 


Small-angle X-ray scattering collimation systems 

ordinarily use slits instead of pinholes in order to give 

sufficient scattered intensity. ¥ith slit collimation, a^ 

a particular angle setting of the collimation system, the 

detector records the scattering from a range of angles, 

rather than from a given angle. This imperfect collimation 

leads to a distortion in the scattering curve. These errors 

are essentially due to the width and height of the direct 

beam. The process of correcting the experimental scattering 

27 - 37 

curve has been considered by several authors . Some of 


the methods suggested need numerical differentiation of the 

experimental curve. Numerical differentiation can be quite 

inexact, specially when applied to an experimental curve 

which is known only approximately. Numerical calculations 

in these methods are tedious and time consuming. A. relatively 

38 

siBiple method has been suggested by Schmidt and Height 
which has the advantage of eliminating the need for numerical 
differentiation or functional fitting of the experimental 
ctarve . For the case of infinite height and negligibly wide 

39 

primary beam this method is readily applicable and is used 
for the present slit-correction of the scattering curve. 
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Small-angle 
ethylene (b) 


X-fay scattering from (a) Poly- 

Irradiated -Poly ethylene (c) ¥ool. 


(d) Standard continuous scale for relative 
intensity measurements* 


PLATE 1 . 




MICROPHOTOMETER RECORD OF LOW- ANGLE X-RAY SCATTERIk 
FROM POLYETHYLENE. 




MICROPHOTOME' 
X-RAY SCATTE 
POLYETHYLENE 




MICRO PHOTOME’ 
X-RAY SCATTE 
POLYETHYLENE 




a 



II, 


STANDARD CONTINEOUS SCALE FOR RELATIVE INTENSITY 
MEASUREMENTS. (MICROPHOTOMETER TRACE) 
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Steps followed in the process are given below. 

For the slits of negligible width and infinite height, 
the experimental scattered intensity F(h) for a given angle 
h is related to the perfect collimation scattered intensity 
l(h) by the relation 


F(h) 


^ OO o 

I{h'^ + 0 d0 


J 


( 5 . 1 ) 


, V . ^ 27,28 

or the scattered intensity l{h) can be written as 

^ _ F’ fch^ + t^)^! (5.2) 

(h^ ^ 


1(h) 


2 

IT 


! 

/ 


L- 


where 0 cind t in equation (5*1) (5.2) are the 

varieties of integration in terms of scattering angle. 

38 

Following the method of Schmidt and Height above 
equation can be expressed as 

OO ~ 

1(h) oc -4-^(2 j+'>)^F(j Ah) - y~ '^ij ^ ((^+j) 


i = 1 


f 5.3'i 


where 


h = 

j Ah 

h. = 

1 

(i+j) A^h 

It 

(i+J) ZiJj 

\ - 

1 

"*i j 

R . = 

1 

Ah (i + 2ia)'^ 


/Vh is taken as 5x1 0"^ radians. Summation for i 
is taken from 1 to 15 which gives fairly good results on a 
test problem. The experimental and slit corrected curves 



for polyethylene anJ irradiated -polyethylene samples are 
shown in Figures 15 and 16 respectively. 

RESULTS 

Studies of loly ethylene Sample 

The slit corrected experimental curve (Figure 15) is 
used for getting the position-correlation function for poly- 
ethylene sample. The relation (4.7) has been used. This 
curve has been plotted and shown in Figure 17* The effect 
of extrapolation of the experimental curve to zero-angle on 
the correlation function has been discusser in the previous 
chapter (p. 66 ), The essential idea is that the relevant 
part of the correlation cxrrve is almost insensitive to the 

extrapolation. Several meixima are obtained. The position 

o . 

of the first maximum at 118 A denotes an approximate repeat 

distance. This distance includes the widths of the crysta- 

ls 

nine and amorphous regions. Buchanann and Miller have 
reported in polyethylene single -lamellar shadow— height s 
varying between 97 and 127 A. The mean shadow-height 
obtained by them from measurements on fracture faces 

o 

several lamellae thick is also found to be nearly 117 A 

which is in general agreement with the result obtained from 

o 

the correlation curve (repeat distance, II 8 A) • 

Correlation-curve can also be used to obtain a very 
rough estimate of the degree of order (0) within the system. 
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1 

This is given by 0 = (l+CD) where CD is the depth of 
the first minimum as shown in Figure 17“ I't gives a value 
of S3fo. This method of estination is based on the assum- 
ption that the distribution of thicknesses in the crystalline 
and amorphous regions is narrow and is peaked around the 
average values. The correlation function is theoretically 
obtained from expression (4.24). The terms upto *^cacac^^) 
are used. The model parameters are the crystallinity (0) , 
the widths of the distribution functions of crystalline and 
amorphous regions ( l\,c and ^a) . The long period (c+a) is 
taken as the position of the first maximum. The parameters 
of the model are varied to match the theoretical curve to 
the experimental correlation function. These curves are 
shown in Figure 1?. This gives a crystallinity of 83^. 

Th© average widths of the crystalline aind amorphous regions 
are 98 and 20 A. The values of the theoretical correlation 
curve with respective terms ( equation (4#24) as a function 

25 

of distance^ are given in Table 7A. Buchanann and Miller 
have reported these distances as 110 and l4 A for 89 ^ degree 
of crystallinity. 

To test the validity of the two-phase model, th© 
irradiated sample of polyethylene has been studied# The 
experimental and slit corrected curve is shown in Figure 16. 
The position-corr elation curve obtained from the slit 
corrected scattering curve is shown in Figure 18* In the 
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TABLE 7 


CALCULATED POSITION -COi®ELAT ION FUNCT ION S* 



C = 

a = 

(a) Polyethylene 

98.0 A Ac = 35.0 

20.0 A Aa = 7.5 

0 

A 

K 

= 83 ^ 


x(A) 

%(-) 

q (x) 

q (x) 

c ac ac ' ’ 


y (x) 

> exp^ '' 

x{A) 

0 

1 .000 

0.000 

0.000 

1 .000 

1 .000 

0 

20 

0.797 

0.030 

0.000 

- 0.021 

0.031 

20 

4 o 

0.599 

0.198 

0.000 

-0.194 

-0.127 

40 

6 o 

0 . 4 i 3 

0.390 

0.002 

-0.154 

-0.155 

60 

80 

0.253 

0.556 

0.009 

- 0.078 

- 0.087 

80 

100 

0.132 

0.672 

0.027 

0.008 

0.012 

100 

120 

0.057 

0.712 

0.070 

0.053 

0.031 

120 

140 

0.020 

0.668 

o.i 47 

0,028 

0.024 

i 4 o 



(B) 

Ir r adia t ed -Po ly e thy 1 ene 



C = 

67.0 A 

/\c = 40.0 

0 

A 







0 ^ 

= 69.9^ 



a = 

29.0 A 

A a = 16.0 

A 



x(A) 

%(*) 

q (x) 

c ac ' ’ 

q (x) 

cacao ^ ’ 

^the^^^ 

' exp^ ' 

x(A) 

0 

1 .000 

0.000 

0-000 

1 .000 

1-000 

0 

20 

0.728 

0.031 

0.000 

0.200 

0-230 

20 

40 

0.486 

0.155 

0.003 

-0.182 

- 0.285 

40 

6o 

0,291 

0.339 

0.014 

- 0.184 

-0.215 

60 

80 

0.152 

0,484 

0.038 

- 0.084 

-0.050 

80 

100 

0.068 

0.546 

0.081 

- 0.014 

-0.060 

100 

120 

0.025 

0.524 

0. l 46 

-0.012 

-0.005 

120 


y, ( x) = Theoretical correlation function; ~ Bxp- 

erimental correlation function. 

All other notations are same as given in text, Chapter IV, 
pages 66-72. 
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present case, the position of the first maximum is not very 

O 

Cleary however, it occurs at rw 84 A. An estimate of the long 
range order is made from the position of the first minimum, 

O 

which gives a distance of So A. Crystallinity obtained 
from the depth of the first minimum 0 = (1 +C'd') ^ gives 
a value of AJ 75^. The theoretical calculations were tried 
by assuming an initial value of crystallinity as obtained 
above and the parameters of the model (0, Ac, A were 
varied. The theoretically calculated curve is shown in 
Figure 18 along with the experimental correlation function. 
The values of the theoretical curve with different terms 
( relation (4.24) ), as a function of distance, are given in 
Table 7B. As expected in an irradiated sample, a low 
crystallinity and wider distribution of widths are obtained. 
The deviation, specially at the first maximum can be inter- 
preted as arising from the changes occurred by irradiation . 
Irradiation brings about cross-linked regions which could 
have a higher electron density and thus two-phase model is 
not quite true. Calculation of correlation function based 
on more phases will yield better results. 

Present analysis is based on the existence of two 
clear-cut phases - the ordered and disordered ones. Such 
a demarcation is not always possible and it is likely that 
a model with continuous variation in electron density would 
yield better results. This and other limitations of the 
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present model are discussed in the concluding chapter of 
the thesis. 


Studies of Wool Sample 

In this part, we briefly describe our preliminary 
results on small-angle scattering from Rambouillet 715 
sample of wool, washed and soxhlet extracted. It was 
obtained from the Western Utilization Research Centre, Albany, 
California. Small-angle X-ray scattering from wool sample 
is given in Plate 1 (c). The microphotometer record of the 
scattering from wool is shown in Figure 19. The analysis is 
carried out by a direct comparison between the scattering 
calculated for certain theoretical models and the scattering 
actually observed in experiment. It is assumed that the 
fibres studied have a high dje.gr.eje of orientation, the micro - 
crystallites being parallel to the fibre axis (nearly so) . 

The scattering obtained consists of a narrow band along the 
equator of about the same width as the height of the beam 
( Plate 1 (c) ). The crystallites may, therefore, be regarded 
as very long compared to their radius and the wavelength of 
the radiation. For the calculation, these particles may be 
considered to have infinite length. The scattering theory 
used is based on a model consisting of a system of patrallel 
cylindrical rods which scatter either independently or 
nonindependently and which have nearly constant diameter 
and homogeneous inner structure. 
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The amplitude 0 scattered by a cylindrical particle 


is proportional to 


2 (kR) 


in which k = , 9 ^ q being half the scattering angle 

and R the radius of the cylinder. For small angles 

^ jy being the scattering angle in radians . 
The function 0(kR) is the form factor of the particle. If 
instead of one cylinder, an assemblage of parallel uniform 
cylinders is present and the elementary waves from single 
particle do not interfere then the total scattered radiation 
is the sum of the radiations from the separate cylinders or 


NjjiCkR)'^ 


in which N is the number of cylinders. The 


square of the form factor has been drawn as a function of 

the angle in radians and matched with the experimental points. 

o 

The curve corresponding to R = 45 A seems to fit the 
experimental curve neairly well. These curves are shown in 
Figure 20. It has been reported that for particles of vary- 
ing size the small-angle scattering is mostly governed by 
the higher end of the scattering periodicity. As the tail 
of the scattering cvcrve does not fit on theoretical curve, 
it is quite possible that there is (a) distribution in the 
size of the particles (b) the interparticle interference 
effects are in^ortant. A crude volume filling calculation 
shows that each microcrystallite unit has approximately 



Intensity ( Arb • Units ) 


0 



* 

FIG. 20 EXPERIMENTAL POINTS AND THEORETICAL 
CURVES' FOR WOOL SAMPLE 



30 alpha helices. An interesting feature of the small- 
angle scattering is that they vary widely from one member 
of the Keratin-myosin group to another even though their 
patterns at wide-angles (i.e. the alpha structure) are 
closely similar. 

Limitations of the present model, specially the effect 
interparticle in"' erf er ence, are discussed in the next chapter. 
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CHAPTER - VI 
CONCLUSION 


The work outlined in the previous chapters is subject 
to certain limitations. For exajuple , vibrations of TATNB 
have been obtained for an isolated molecule though in practice 
because of the insolubility of TATNB in almost all the solvents, 
it is just not possible to obtain the spectra in solution form. 
The experimental data reported here is, thus, obtained in 
solid form. However, physical state of affairs in a crysta- 
lline form is different from the sinqslified model of isolated 
molecule . A complete interpretation of the vibrational modes 
would require calculations on a three dimensional system, 
taking into account all kinds of intermolecular interactions. 
The size of the dynamical matrix involved in the solution of 
vibrational problem of a three-dimensional system is very large 
and makes the solution of secular equation almost prohibitive. 
The general effects on the frequency spectrum to be 

expected as a result of intermolecular interactions in a 

1 

crystal can be fairly easily visualised . According to 
2 

Halford the vibrations of the molecule in the crystalline 
state are governed by new selection rules derived from site 
symmetry - a local symmetry around the centre of gravity of 
a molecule in a unit cell. In general, site symmetry is 
lower than the molecular symmetry in am isolated state. The 
change in symmetry may split the degenerate vibrations and 
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activate optically inactive vibrations. In addition, the 
spectra obtained in the crystalline state exhibit lattice 
modes. A more complete analysis including lattice modes 
can be made by the factor group analysis corresponding to 
the crystal space group. For a unit cell containing N 
atoms there will be (3N-3) normal vibrations. For the 
case of TATNB , where there are two molecules in a unit cell, 
total number of normal vibrations will be l4l (2x24x3 - 3) * 
The 12 fundamentals known as external vibrations can be 
further divided into nine rotatory type and three transla- 
tory type of oscillations. A complete analysis of factor 
group fundamentals, however, can be carried out only when 
the nature of the intermolecular interactions is known. 

Assuming a admittedly simplified model of isolated 
molecule, it has, nevertheless, been possible to find 
reasonable correlation with observed spectra. Although, 
interactions through hydro gen -bonds have not been explicitly 
considered, the adjustment of force constants to obtain 
agreement with observed frequencies has effectively ta en 
into account these interactions in a circuitous manner. 

The force field obtained cam be used for similar molecules 
as a starting point, since the validity of transferability 

of Urey -Bradley force field has been shown in a series of 

3,4 


papers by Scherer 
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Coming to the second part of the thesis, i.e*, the 
work on small -angle scattering of X-rays, the two-pfaase 
parallel layer model assumed for polyethylene sample seems 
to hold to a fair degree. This is seen from the fitting of 
the major part of the experimental correlation curve with 
that of the theoretical one. However, small deviations are 
observed at the first minima. Similar deviations are observed 
in case of irradiated polyethylene sample where the first 
maximum is neither well defined nor well matched with the 
theoretical curve. The situation in irradiated sample is 
different because irradiation brings about cro ss -linking^ * 

The analysis of the experimental curve is complicated 
because the cross-linked regions could have a different 
average electron density than the crystalline and amorphous 
regions. This will introduce a coirpletely different phase* 
Irradiation also brings about the broadening of the distri- 
bution of layer thicknesses. This suggests that for an 
irradiated sample a model with more than two phases is 
needed ♦ 

The present analysis of the experimental curve for 
polyethylene is based on the existence of two cleair-cut 
phases - the ordered and the disordered ones. Such a 
demarcation is generally not possible. It is likely that a 
model with continuous variation in electron density at the 
phase boundaries would yield better results. The correlation 
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function approach for melt-crystallized polymers has several 

advantages as indicated in previous chapters. The success 

of the two-phase model suggests that it can be used for 

similar systems. Recently, the molecular parameters, such 

as average length of the crystalline and amorphous phases 

have been obtained for semicrystalline polymers by making 

5 7 

use of the mathematical methods of Markoff chains ’ . It 

would be interesting to compare the present work with such 
calculations and obtain transition probability of various 
states. It should be possible to relate these probabilities 
with the nature of potential barrier around single bonds* 

A preliminary analysis of the scattering curve from 
wool sample is reported in the thesis. The conqplete analysis 
has to include interparticle interference which plays an 
important role in these systems. From the experimental 
point of view these effects can be removed by swelling the 
sample in a suitable swelling agent. Further work, taking 
into account the above factors, is very essential for a 
fuller interpretation of the experimental results. 
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APPENDIX 


SYMMEIRY COORDINATES 
FOR TATNB ^ 

Definition of Internal Coordinates 


Stretches 


to 

^6 


r(C-C) 


"7 

to 

^12 

== 

r (N— 0 ) 


^13 

to 

^18 

= 

r(N-H) 


^19 

to 

^21 

=: 

r{C-S)Q 


^22 

to 



r(C-N)g 

In-plane angle bendings 

^25 

to 

^30 

= 

0(C-C-«)q 


^31 

to 

^36 

= 

0 ( CHSF— 0 ) 


^37 

to 

^42 


0(C-C-N)jj 



to 

«48 

= 

0(CHN-H) 



to 


= 

0(C-C-C)q 


^52 

to 

^54 

= 

0( 0-^-0) 


R55 

to 

R37 

= 

0(C-O-€)jj 


^58 

to 

^60 


0(H-^H) 

Torsional coordinates 

^61 

to 

^66 

= 

t(C-C) 


^67 

to 

^69 

= 

t(C-N)Q 


^0 

to 

®72 


t(C-N)jj 

Out -of -plane wagging 

^73 

to 

R75 

= 

tjCOgN-D) 


^76 

to 

^78 

= 

^^(HgN-C) 


^79 

to 

^81 

= 

10{B -N)q 


^82 

to 

^84 

= 

U(B-N)jj 


* 'B' represents benzene ring. 

r stands for bond stretching, 0 for in-p\ane angle 
bending, t for torsions and for out-of -plane vags . 



In-plane symmetry coordinates 


(1) 

vy^ (R^ 

. R2 

(2) 

1//^ (R^ 

+ Rg 

( 3 ) 

1/76 (R^3 

+ ^14 

(M 

V /3 (R^9 


( 5 ) 

vys (R22 

■*■ ^23 

(6) 

1//^ (R23 

^26 

( 7 ) 

v/^ (R3, 

^32 

(8) 

1/-/6 (R3^ 

+ R3Q 

( 9 ) 

1//6 (R43 

+ 1^44 

(10) 

1/73 (R^3 

^50 

(n) 

1/73 (R32 

+ 

(.IS) 

1/73 (R33 


( 13 ) 

1//3 (R38 

+ R59 


(1) 

1/7^ 


-^2 

(2) 

1/76 

(■'t 


(3) 

1/76 

<‘'13 

- ^14 

(^) 

1/76 

(1*25 

' ^26 

(5) 

1/76 


“ ^32 

(6) 

1/7^ 

(R3, 

“ ^38 

(7) 

1/7^ 

(«43 

~ ^44 


Species 

+ + R^ + ) 

^9 ^ ^10 ®11 ^ 12 ) 

+ R^^ + + R^g) 

+ ^^ 24 ) 

+ R^^ + R 2 g + R 29 + R 3 Q) 
+ R^^ + R^^ + R^^ + R^^) 

+ ^ 3 ^ + R^^q + R^^ + R^ 2 ) 

+ Rj^^ + R^^ + R^^ + R^g) 

.R 3 ,) 

+ R54) 

+ « 6 o) 


Species 


+ 

-R4 

+ 

>*5 

- «6 ) 


" ^10 

4 - 

«11 

- "12) 

♦ R,, 

- ^16 

+ 

®17 

' ® 18 ^ 

+ R2y 

' ^28 


^29 

' ® 30 ^ 

+ R33 

- R34 

+ 

^35 

- 

+ R39 

“^0 

+ 

^41 

- ^42) 

+ ^45 

■ ^46 

+ 

R47 

- ^ 48 ^ 
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Species 


(T) 

1//12 

(2R^ 

♦ ®2 

(2) 

I//T2 

(2R^ 

- "'*2 

( 3 ) 

1/2 

( R3 

+ R^ 

(^) 

1/2 

( R3 

- «4 

( 5 ) 

i/yT 2 

(2R^ 

. 2Rg 

( 6 ) 

1//T2 

(2R^ 

- 2^8 

( 7 ) 

1/2 

( H, 

"■ ^10 

( 8 ) 

1/2 

(Kj 

-^10 

( 9 ) 

i/yT2 

( 2 R ,3 

+ 2R^^ 

( 10 ) 

i/yl 2 

( 2 E ,3 

- 2^14 

( 11 ) 

1/2 

( R ,5 

^ ^16 

( 12 ) 

1/2 

( Rij 

- ^16 

(13) 

V/^ 

( 2 R ,5 

^20 

( 14 ) 

1//2 

( «20 

^21 ) 

(^ 3 ) 

1//6 

(=«22 

■ ^23 

(16) 

1//2 

( ^23 

- « 24 ) 

(17) 

1//‘T2 

( 2^25 

- "^26 

( 18 ) 

1//T2 

( 21*25 

- ^^26 

(19) 

1/2 

{ R27 

^28 

(20) 

1 /2 

( *<27 

" ^28 


_ Rg ) 

R3 + Rj^ - R3 + Rg ) 

^5 - ""6 ) 

S ""6 ) 

^ -""lo -^11 -^12) 

^9 ^ ^10 ■ ^11 + ^12) 

"'ll - ^12) 

^11 " ^12) 

^15 " ^^16 ■ ^17 - ^is) 
^15 ^16 " ^17 + ^ 18 ^ 

^17 ‘ ^18^ 

^17 ®18^ 

Rgi ) 

^24) 


^27 

t 

00 

\D 

0 

fE| 

1 

^27 

^28 " ^29 

^30 

^29 

1 

pa 

0 


^29 

+ R30) 




(21) 

1//12 

( ^ + 

2R32 - R33 - R34 - R35 - ^35) 

(22) 


(2R3, - 

2R32 - R33 + f ^34 - ^35 + ^36^ 

( 23 ) 

1/2 

( R33 * 

^34 ■ ^35 “ ^36^ 

( 24 ) 

1/2 

( R33 - 

^34 " *^35 ^36^ 

( 25 ) 

i/yT2 


2R38 ~ ^39 " ^ 4 o ” ^ 4 i ” ^ 42 ^ 

(26) 

I//T2 

(2R^^ - 

2R38 - ^39 + ^ 4 o ■ ^41 ^42) 

( 27 ) 

1/2 

( K39 t 

^ 4 o “ ^ 4 i “ ^ 42 ^ 

( 28 ) 

1/2 

( «39 - 

^ 4 o " ^41 ■*■ ^ 42 ^ 

( 29 ) 

I//T2 

( 2 R ^3 + 

2R44 - R43 - ^45 “ ^47 ■ ^ 48 ^ 

( 30 ) 

i/yT2 

(®43 - 

2R44 - R45 ^^46 - ®47 ^ ^ 48 > 

( 31 ) 

1/2 

( R43 + 

^46 “ ^47 *" ^ 48 ^ 

( 32 ) 

1/2 

( R45 - 

®46 ” ^47 ® 48 ^ 

( 33 ) 

1//6 

(2H49 - 

^50 “ ^51^ 

( 34 ) 

1//^ 

^ ^^50 " 

R51) 

( 35 ) 

1//^ 

(2R52 - 

R53 - R54) 

( 36 ) 

1//2 

( R53 - 

R54) 

( 37 ) 

1//6 

(2R55 - 

R56 - ^ 57 ^ 

( 38 ) 

i/ys 

( «56 - 

R57) 

( 39 ) 

1//6 

(2R58 - 

^59 - ^60) 

( 40 ) 

1/72 

( 1^59 - 

^^6o) 





Out -of -plane symmetry coordinates 




If 

Specxes 

(1) 

■\//Z 

( "61 * 

^62 ®63 ^64 ^65 ^66^ 

(2) 

1//3 

( "67 * 

^68 ^69^ 

( 3 ) 

1//3 

{ 

^71 ^72^ 




Species 

(1) 

1/y^ 

( "61 - 

®62 + ^63 " ^64 + ^65 ^ ^66^ 

(2) 

i/yi 

( ^>73 * 


( 3 ) 

i/y 3 

( + 

4 - R,^g) 

(^) 

i/yi 

( ^79 + 

^80 ^ 81 ^ 

( 5 ) 

1/73 

( Rg2 + 

^83 ^ 84 ^ 




e" Species 

(1) 

1//T2 

(2^61 ^ 

2^62 - ^63 - ""64 - ^65 " ^ 66 ^ 

(2) 

1/A2 

(21^61 - 

2R^2 - ^^63 ^ ^64 - ^65 ^ 66 ) 

(?) 

1 /2 

( R 63 + 

^64 “ ^65 " ®66) 

(^) 

1/2 

( ^63 - 

^64 ■ ^65 ®66^ 

( 5 ) 

l/v/6 

(^^61 - 

^68 “ ^69^ 

(6) 

1//2 

( ^68 - 

^69) 

( 7 ) 

1//6 

(2R^0 - 

^71 “ ^72^ 

(8) 

1//2 

( Ry^ - 

R72) 

( 9 ) 

1//^ 

( 2 Ry 3 " 

R74 - J^ 75 > 

(10) 

1//2 

( ^74 

R75) 

(11) 

1//6 

(^^76 

Rijrjr ~ ^78^ 

(12) 

1//2 

( 1^77 - 

1^73) 

( 13 ) 

1//6 

( 21^79 “ 

®80 - "" 81 ) 

( 14 ) 

lX/2 

^ ^80 

^81 ) 

( 15 ) 

1/y^ 

( 2 R 82 - 

^83 " ^ 84 ^ 

(16) 

1//^ 

( ^^83 ■■ 

^ 84 ) 
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